From constructive field theory to fractional 
stochastic calculus. (I) The Levy area of fractional 

Brownian motion with Hurst index a 6 (hi) 



Jacques Magnen and Jeremie Unterberger 



Let B = (Bi(t), . . . ,Bd{t)) be a d-dimensional fractional Brownian motion 
with Hurst index a < 1/4. Defining properly iterated integrals of B is a 
difficult task because of the low Holder regularity index of its paths. Yet 
rough path theory shows it is the key to the construction of a stochastic 
calculus with respect to B, or to solving differential equations driven by B. 

We show in this paper how to obtain second-order iterated integrals as the 
limit when the ultra-violet cut-off goes to infinity of iterated integrals of 
weakly interacting fields defined using the tools of constructive field theory, in 
particular, cluster expansion and renormalization. The construction extends 
to a large class of Gaussian fields with the same short-distance behaviour, 
called multi-scale Gaussian fields. Previous constructions [361 135] were of 
algebraic nature and did not provide such a limiting procedure. 
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Introduction 

A major achievement of the probabilistic school since the middle of the 20th 
century is the study of diffusion equations, in connection with Brownian mo- 
tion or more generally Markov processes - and also with partial differential 
equations, through the Feynman-Kac formula - with many applications in 
physics and chemistry |38j . One of the main tools is stochastic calculus with 
respect to semi-martingales M . An adapted integral such as X(u)dM(u) 
may be understood as a limit in some sense to be defined. Classically one 
uses piecewise linear interpolations, X>s<t 1 <...<t JV <t 

X(ti)(M(t i+1 ) - M(U)) 

or ^2 s <t!< <t N <t X + * (M(ti+i) — M(ti)); these approximations de- 
fine in the limit N — > oo the ltd, resp. Stratonovich integral. The latter one 
is actually obtained e.g. if M = W is Brownian motion and X(t) = f(Wt) 
with / smooth as the limit lim £ _>.o J f(W s (u))dW £ (u) for any smooth ap- 
proximation (W e ) e >o of W converging a.s. to W (see [3D] , or [T7] p. 169). 
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The Stratonovich integral J s * X(u)d Strato M(u) has an advantage over the Ito 
integral in that it agrees with the fundamental theorem of calculus, namely, 
F(M(t)) = F{M(s)) + JlF'{M(u))d strato M{u). 

The semi-martingale approach fails altogether when considering stochas- 
tic processes with lower regularity. Brownian motion, and more generally 
semi-martingales (up to time reparametrization), are (l/2)~-H61der, i.e. a- 
Holder for any a < 1/2 0. Processes with a-H61der paths, where a <C 1/2, 
are maybe less common in nature but still deserve interest. Among these, 
the family of multifractional Gaussian processes is perhaps the most widely 
studied [28], but one may also cite diffusions on fractals [15], sub- or su- 
perdiffusions in porous media [11] [18] and the fascinating multi-fractal ran- 
dom measures/walks in connection with turbulence and two-dimensional 
quantum gravity [H |6] . Although some of these processes appear in applica- 
tions rather as a random landscape in D > 1 space dimensions (which raises 
related but different problems), we view them here - restricting ourselves 
to D = 1 - as some colored noise integrated in time, and wish to define 
stochastic integration with respect to them or to solve differential equations 
driven by them. 

We concentrate in this article on multiscale Gaussian processes (the 
terminology is ours) with scaling dimension or more or less equivalently 
Holder regularity a G (0,1/2), the best-known example of which being 
fractional Brownian motion (fBm for short) with Hurst index a, B a {t) 
or simply B(t) q We consider more precisely a two-dimensional fBm, 
B(t) = (Bi(t), B 2 (t)), with independent, identically distributed compo- 
nents El- The simplest non-trivial stochastic integral is then A(s,t) := 
J* s * dBi(t\) f^ 1 dB 2 (t 2 ) = f^(B 2 (u) — B2(s))dBi(u), a twice iterated integral, 
where B = (B\(t), B^t)) is a two-component fBm. Since f* dBi(ti) J"* 1 ^-62(^2)+ 
r s dB 2 {t 2 ) j'l 2 dB^tx) = {B 1 {t)-B 1 {s)){B 2 (t)-B 2 (s)), one is mainly inter- 
ested in the antisymmetrized quantity (measuring an area, as follows from 

1 Recall that a continuous path X : [0, T\ — > R is a-Holder, a £ (0,1), if 

SUP s , te[0 ,T] | t - 3 |c Sl < OO. 

2 It is (up to a constant) the unique self-similar Gaussian process with stationary in- 
crements. The last property implies that its derivative is a stationary field. 

3 The one-dimensional case is very different and much simpler, and has been treated in 

[El- 
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the Green- Riemann formula), 

CA(s,t) := f dfli(ti) I' dB 2 (t 2 ) - f dB 2 (t 2 ) f ' dBifa) 

= J\b 2 (u) - B 2 («))dBi(«) - (Sx(«) - B x (s))dB 2 (u), 

(0.1) 

called Levy area. The corresponding Stratonovich integral, obtained as a 
limit either by linear interpolation or by more refined Gaussian approxima- 
tions [5j [271 Ell [33], has been shown to diverge as soon as a < 1/4. 

This seemingly no-go theorem, although clear and derived by straight- 
forward computations that we reproduce in short in section 1, appears to 
be a puzzle when put in front of the results of rough path theory [231 [241 [T3"l 
[20]. [2l~] 110] . We shall not enter into the details of this fascinating theory, 
with its deep connections with nilpotent groups, Hopf algebras and sub- 
Riemannian geometry, leaving this to a forthcoming paper |25] where these 
structures will be instrumental to define higher-order integrals. Let us only 
state that rough path theory allows one to define integration with respect 
to an a-H61der path T in terms of substitutes of iterated integrals of T up 
to order |l/ a J) L • J=entire part, called rough path over T, satisfying some 
Holder regularity property and also two algebraic properties, called respec- 
tively Chen and shuffle property. The former one is of geometric nature and 
may be seen as an additivity property for the areas or volumes generated 
by double or multiple iterated integrals. In particular, it implies the follow- 
ing compatibility property between the twice iterated integral A and the 
original process B, 

A(s, t) = A(s, u) + A(u, t) + (B 2 (u) - B 2 (s))(B 1 (t) - B x (u)). (0.2) 

The shuffle property, on the other hand, ensures once again that the funda- 
mental theorem of calculus holds. Whatever the exact formulation of these 
two properties, they hold true trivially for the usual iterated integrals of T 
if r is smooth, and (being of algebraic nature), also hold true trivially when 
the rough path over T is constructed out of true iterated integrals of smooth 
paths by some limiting procedure as explained above. 

The main point is now the following. General theorems show (i) that 
rough paths always exist (existence theorem); (ii) that any roughpath over 
r may be constructed by some rather abstract limiting procedure □ (approx- 
imation theorem). Unfortunately, these theorems are not constructive; also, 

4 using pieces of horizontal sub-Riemannian geodesies. 
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they are of pathwise nature, and thus not necessarily appropriate in the 
probabilistic setting which is our starting point. Recent results using the 
underlying Hopf algebra structures [361 E3 EH El EI] give an explicit, gen- 
eral solution to the existence problem, well-suited in particular for Gaussian 
processes, but no one knows how to obtain these rough paths by an explicit 
limiting procedure. 

Our project in this series of papers is to define an explicit rough path over 
fBm with arbitrary Hurst index, or more generally multiscale Gaussian fields 
(to be defined below) by an explicit, probabilistically meaningful limiting 
procedure, thus solving at last the problem of constructing a full-fledged, 
Stratonovich-like integration with respect to fBm. 

Roughly speaking, it is obtained by making B = (B{\),B(2)) interact 
through a weak but singular quartic, non-local interaction, which plays the 
role of a squared kinetic momentum, associated to the rotation of the path, 
and makes its Levy area finite. Following the common use of quantum field 
theory, this is implemented by multiplying the Gaussian measure by the 
exponential weight e -§«4//An*(*i,*a)(ti,*a)|ti-te|- 4a cfti* I g j with 

Ar^i, 02) (*i,* 2 ) = A 2 {(dA + )(h)(dA + )(t 2 ) + (dA-)(h)(dA-)(t 2 )} , 

(0.3) 

where: A (the coupling parameter) is a small, positive constant; <pi,^>2 are 
the (infra-red divergent) stationary fields associated to B\,B 2 , with covari- 
ance kernel as in eq. (|2.2ip . and similarly, A ^ are stationary left- and 
right-turning fields, built out of (pi, (p 2 , and representing the singular part 
of the Levy area (see section 1 for details). By assumption a < 1/4, so that 
the kernel \ti — t 2 \~ 4a is locally integrable. The statistical weight is maximal 
when dA + = dA~~ = 0, i.e. for sample paths which are "essentially" straight 
lines. Another way to motivate this interaction is to understand that the 
divergence of the Levy area is due to the accumulation in a small region 
of space of small loops |20j ; the statistical weight is unfavorable to such an 
accumulation. On the other hand, the quantities in the first-order Gaus- 
sian chaos (in other words, the n-point functions (Bi 1 (xi) . . . Bi n (x n ))\ = 



IE 



B h (xi) . . . Bi„ (x n )e-K // £inMxM){txte)\tx-H 4a d tl dt 2 



where Z := E 



,-h C 'cf I £int(<l>l,(l>2){tl,t2)\tl-t2\- 4a dt 1 dt2 



^1) • • • > — l)^) 



is a normalization con- 



stant called partition function) - associated to pure translation - are insen- 
sitive to the interaction. Thus the compatibility property (|0.2p is satisfied, 



5 The unessential constant c' a is fixed e.g. by demanding that the Fourier transform of 
the kernel c4|ti - t 2 \~ ia is the function |f| 4a_1 . 



5 



and one has constructed a Levy area for fBm. 

It would be very interesting to find an equivalent pathwise description 
of this interaction in terms of a constrained, long-memory self-interacting 
walk, see section 1 for more comments. The interaction Ci n t would thus 
play the role of a Ginzburg-Landau functional, whose usual, local version 
is conjectured to give e.g. a field-theoretic description of the Ising model. 
Proving such an equivalence should be much easier here (despite the unusual 
non-local features) since it is a one-dimensional model. 

Starting from the above field-theoretic description, the proof of finiteness 
and Holder regularity of the Levy area for A > small enough follows the 
well-established scheme of constructive field theory. The main theorem may 
be stated as follows. As a rule, we denote in this article by E[...] the Gaussian 
expectation and by {...)\ tP the expectation with respect to the A-weighted 
interaction measure with scale p ultraviolet cut-off, so that in particular 

E[...] = (...) 0l OQ. 

Theorem 0.1 Assume A G (^,\)- Consider for A > small enough the 
family of probability measures ( also called: ((f), d(f>, cx)-model ) 

P^icpufo) = e-5 c '-// d *i d *2l*i-*2l~ 4Q ^(*i^)(*i'*2)d^(0 1 )d/i^( ( />2), 

(0-4) 

where dp~^ p (<j)i) = dp,(<j) i p ) is a Gaussian measure obtained by an ultra- 
violet cut-off at Fourier momentum |£| ~ M p (M > 1), see Definition \2.10\ 
and section 5. Then (P^ p ) p converges in law when p — > oo to some mea- 
sure P\, and the associated twice iterated integral f d(j)^ p (ti) J** 1 d(f)^ p (t2) 
converges to a well-defined quantity A(s, t) satisfying the Chen and shuffle 
properties, with variance 

(A(s, t) 2 ) x := lim (A(s, t) 2 ) A , p = (i^l + K 2 )\t - s\ 4a . (0.5) 

p^oo A z 

Furthermore, the connected moments of A{s,t) in the interacting theory 
- except its variance - coincide with those of the free, Gaussian theory, allow- 
ing for the computation of moments of A(s,t) in terms of explicit Gaussian 
integrals. In particular, one finds: 

{A(s, t) 2n ) x < C n \t - s\ 4na , n > 0. (0.6) 

As an immediate corollary following from the Kolmogorov- Cents ov lemma 
A(s,t) has a.s. (2a) -Holder paths. 

6 This lemma states that a process X such that E\X(t) - X(s)\ p < C\t - s| 1+p/3 has 
almost surely /3~-H61der paths. Classical hypercontractivity arguments implies the same 
conclusion for a Gaussian process provided simply E\X(t) - X(s)\ 2 < C\t - s\ 2 > B . 
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The only difficulty is to prove the convergence of the moments of finite 
distributions of the Levy area when p — > oo; let us see how to deduce the 
assertions of the theorem from this. 

The weak convergence of (some subsequence of) the family of laws (P7 >p )p 
is actually a consequence - using standard arguments including Prohorov's 
theorem, [2], chap. 6, 8, 12 - of the following estimates which imply that 
(P^")p is tight, 

<|0(t)-0( S )| 2 ) A , p <#|f- S |^ (0.7) 

K being a constant uniform in p. This last estimate is proved in a stronger 
Fourier form0, 

(1^)1 V < T^tW (0-8) 

Finally, standard multi-scale estimates, sketched in section 1, show that 
the generating function of connected moments of finite distributions of the 
Levy area, («i, . . . ,u n ) h-> log(exp i (u\A{si, t\) + . . . + u n A(s n , t n )))\ )00 , has 
a non-zero radius of convergence, hence (see [7], Chap. XV, §4) the mo- 
ments of the finite distributions of the Levy area determine its law, and 
every subsequence of (P^ p ) p converges to the same limit. It would also be 
possible (by a simple extension of the results of section 6) to prove a speed 
of convergence for these moments, e.g. (s)<^ (t)} \ ;P — {<j>> (s)^ (t))\ iP > = 
Kr ? 3 2 \t"l s | y O ( M ~ ^ 8a ~ 1 ~^ ) (f° r some constant K r depending on r = 
1,2,...) if j < p' < p. 

The result is not difficult to understand using the non-rigorous perturba- 
tion theory. First, by a trick explained in section 1, one replaces the non-local 
interaction £(0i, 4>2)(ti, ^2)^1 — ^2| _4a with a local interaction C(<p\, 4>2,cr)(t) 
depending on a two-component exchange particle field a = (<r + (i), o"_(i)). 
Then a Schwinger-Dyson identity (a functional integration by parts) relates 
the moments of A to those of a. Simple power-counting arguments show 
that a connected 2re-point function of a alone is superficially divergent if 
and only if 1 — Ana > 0. Thus, restricting to a > 1/8, one only needs to 
renormalize the two-point function. Since the renormalized propagator of a 
is screened by an infinite mass term (as we shall see in section 1), the theory 
is free once one has integrated out the <r-field, hence one retrieves the under- 
lying Gaussian theory (</>i,</>2). The Schwinger-Dyson identity then shows 
that the two-point functions of A have been made finite. 

7 Namely, |0(t) - 0( S )| 2 < /, e|> ^ + k\<j^ K %f+$ d *> < K '\ l ~ s | 2 ° '• 
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Since constructive field theory \26\ [30] is not a standard tool in prob- 
ability, we shall spend some time explaining in full details and generality 
its objects, in particular the cluster expansions which lie at its heart. Clus- 
ter expansions are based on a simplified wavelet decomposition (V'a) °f the 
fields V = 01 ? 02 ) °"± where j is a vertical (Fourier) scale index, and A a 
horizontal (i.e. in direct space) interval of size M~ 3 around the center of 
the wavelet component. Each ifx^ is to be seen as a degree of freedom of 
the theory, relatively independent from the others, so that the interaction 
may be expressed as a divergent doubly-infinite vertical and horizontal sum. 
The horizontal (H) and vertical (V) cluster expansions allow to rewrite the 
partition function Zy P over a finite volume, with ultraviolet truncation at 
scale p, as a sum, 

V = E^ E F HV (F 1 )...F HV (T n ), (0.9) 

n Pi,...,P n non— overlapping 

where: 

- Pi,...,P n are disjoint polymers, i.e. sets of intervals A connected by 
vertical and horizontal links; during the course of the expansion, the Gaus- 
sian measure has been modified so that the field components belonging to 
different polymers have become independent; 

- Fffv(P), P = Pi, ...,P n is the A-weighted expectation value, F HV (F) = 
(fHvi^^x, of some function fnv depending only on the field components 
located in the support of P. 

The fundamental idea is that (i) the polymer evaluation function F(¥) 
is all the smaller as the polymer P is large, due to the polynomial decrease 
of correlation at large distances (for the horizontal direction), and to power- 
counting arguments developed in section 4 for the vertical direction, leading 
to the image of horizontal islands maintained together by vertical springs; 
(ii) the horizontal and vertical links in P (once one interval belonging to P 
has been fixed) suppress the invariance by translation, which normally leads 
to a divergence when \V\ — )■ oo. A classical combinatorial trick, called Mayer 
expansion, allows one to rewrite eq. (|0.9p as a similar sum over trees of poly- 
mers, also called Mayer- extended polymers and denoted by the same letter P, 
but without non-overlap conditions, Zy p = Y^ n £p a p„ ■ ■ ■ F{¥ n ), 

where F = Fhvm is the Mayer- extended polymer evaluation function, so 
that lnZy p = £ P F(P). In the process, local parts of diverging graphs have 
been resummed into an exponential, leading to a counterterm in the inter- 
action; this is the essence of renormalization. On the whole, one finds that 
in the limit \V\, p — > oo, the free energy In Zy p is a sum over each scale 



8 



of scale-dependent extensive quantities, i.e. In Zy p = \V\ Yljtfo ^fv^ P i 
where fy* p converges when \V\ — > oo to a finite quantity of order 0(A). 
One retrieves the idea that each interval of scale j contains one degree 
of freedom. Finally, ra-point functions are computed as derivatives of an 
external-field dependent version of the free energy. 

Here is the outline of the article. Section 1 recalls the classical results 
on the divergence of the Levy area for a < 1/4 and recasts them into nota- 
tions which are more appropriate for cluster expansions. It also contains a 
heuristic section using perturbative field theory and giving some intuition on 
why adding the interaction term d n t briefly described above should make 
the Levy area convergent. Section 2 contains the definition of multiscale 
Gaussian fields, including fBm, and many useful notations and general re- 
sults concerning the scale (wavelet) decompositions. Section 3 is on cluster 
expansions, section 4 on renormalization. Sections 2 to 4 are extremely gen- 
eral, valid also for fields living on R-°, in the hope that they may serve as 
a basis for future work, possibly also for D > 1-models. Section 5, on the 
contrary, concentrates on the definition of our specific ((f), d(j), <r)-model. The 
proof of finiteness of n-point functions of the Levy area and freeness of the 
0-ficld is given in section 6. It depends on Gaussian bounds which hold in 
great generality, and on domination arguments which are very specific of our 
model. 

Notations. Cluster expansions imply the use of many indices and let- 
ters. Let us summarize here some of our most important conventions, in the 
hope that this will help the reader not to get lost. 

1. Quite generally, ip = (ipi(x), . . . ,-0d( x )) * s a d-dimensional field living 
on R D , D > 1, with scaling dimensions . . . , (3d- Roughly speak- 
ing, for probabilists, (3 is the Holder continuity index, at least when 
/3 > 0; physicists usually call scaling dimension —(3. Most of the 
model-independent results here are valid for arbitrary D. The Fourier 
transform is denoted by T . 

2. If tp = ip(x) is a Gaussian field, then its scale decomposition (see 
section 1) reads ip = J2j>o' l P' 1 - The low-momentum, resp. high- 
momentum field with respect to scale j is denoted by ip~^i = Ylh<j 
resp. ipi^ = Xlfc>j ^ ' ■ All through the text, we observe the following 
convention: if h,j,k are scale indices, then h < j < k; any primed 
scale index (for instance j' , p' , . . .) is less than the original scale index 



9 



j, p, ■ ■ ■ Secondary fields (i.e. low-momentum fields tp minus their aver- 
age) are denoted by Sip. We also introduce restricted high-momentum 
fields, see section 1, denoted by Res ip. 

3. (products of fields) If tp = (ipi(x), . . . ,tpd(x)) is a ci-dimensional field, 
and / = (i±, . . . , i n ) £ {1, . . . , d} n (n > 1) is a multi-index, we denote 
by tpi the product of fields ipi{x) := ip^ (x) . . . ipi n (x). The interac- 
tion Lagrangian Ci n t is written in general as YT q =i Kq\ Kq, 4)i q for some 
constants K q and exponents K q . Cluster expansions produce propaga- 
tors and products of fields which are linear combinations of monomials 
(also called: G-monomials) , generically denoted by G. 

4. (constants) K is a constant depending only (possibly) on the details of 
the model, such as the degree of the interaction, the scaling dimension 
of the fields... It may vary from line to line. M > 1 is the basis of the 
scale decomposition; it is an absolute constant, whose value is unim- 
portant. 7 is some constant > 1, or sometimes simply > 1. N ex t t max is 
such that every Feynman diagram with > N extjTnax external legs is su- 
perficially convergent; this perturbative notion also plays a central role 
in constructive field theory. Our model (<p, d<p, a) has N ext:Tnax = 4. 

5. (variables) The maximum scale index is p. Other scale indices are 
denoted by h, j, k, or any of those with primes. Indices i are summation 
indices, with finite range, used in various contexts, for instance for the 
field components. The parameters of the horizontal, resp. Mayer, 
resp. vertical (also called momentum-decoupling) cluster expansion, 
are denoted by s, resp. S, resp. t. The scale of an interval A is denoted 
by j(A). If A is an interval, then n(A) is the coordination number 
of A inside the tree defined by the horizontal cluster expansion, while 
N(A), resp. iVj(A) is the total number of fields, resp. the number of 
fields ipi located in the interval A. r stands for a number of derivations, 
usually with respect to the t-parameters, in which case it is assumed 
to be < N ext)inax + <3(n(A)). 

1 Statement of the problem and heuristics 

The quantity we want to define in the case of fractional Brownian motion is 
the following. 

Definition 1.1 (Levy area) The Levy area of a two-dimensional path T : 
M. — t- IR. 2 between s and t is the area between the straight line connecting 
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(ri(s),r 2 (s)) to (ri(t),r 2 (t)) and the graph {(ri(u),r 2 («));a < u < t}. It 
is given by the following antisymmetric quantity, 



CAr(s,t) := / dTi(ti) / dT 2 {t 2 ) - dT 2 {t 2 ) dT^h). (1.1) 
1.1 A Fourier analysis of the Levy area 

In order to understand the analytic properties of the Levy area of fBm, we 
shall resort to a Fourier transform. One obtains, using the harmonizable 
representation of fBm introduced in subsection 2.2, 

(1-2) 

The Levy area CA(s, t) := CAb(s, t) is obtained from this twice iterated 
integral by antisymmetrization. Note that CA(s, t) is homogeneous of degree 
2a in \t — s\ since B(ct) — B(cs), c > has same law as c a (B(t) — B(s)) by 
self-similarity. 

Expanding the right-hand side yields an expression which is not homo- 
geneous in £. Hence it is preferable to define instead the skeleton integral, 
depending only on one variable, 



i r dWi 
2^c~J \i^- 1 ' 2 U a - 1 ' 2 [i(£i + &)][&]' 



dVFi(6)^ 2 (e 



2 J 



(1.3) 



where by definition J e iU ^du = From „4(i) and the one-dimensional 
skeleton integral 

Mt) = (^c a )-z J dB t (u) = J ■ — , (1.4) 

which is the infra-red divergent stationary process associated to B, see sec- 
tion 2, one easily retrieves A(s, t) since 

A(s,t) = £ dB^h) dB 2 (t 2 ) - J* dB 2 (t 2 )J 

= A{t)-A(s)+A d (s,t), (1.5) 

where (2irc a ) zAd(s, t) := (B\(t) — B\(s))(f) 2 (s) (called boundary term) is a 
product of first- order integrals. 
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One may easily estimate these quantities in each sector |6| ^ |6|- in 
practice, it turns out that estimates are easiest to get after a permutation 
of the integrals (applying Fubini's theorem) such that (for twice or multi- 
ple iterated integrals equally well) innermost (or rightmost) integrals bear 
highest Fourier frequencies; this is the essence of Fourier normal ordering 
[361 E E3- This gives a somewhat different decomposition with respect to 
CDS since f* cLBi(ii) /** dB 2 {t 2 ) is rewritten as - f* dB 2 (t 2 ) // 2 d-Bi(ii) in 
the "negative" sector |6| > \£ 2 \. After some elementary computations, one 
gets the following. 

Lemma 1.2 Let 

A + {t) := 2Trc a j dt x J'' dt 2 F- 1 ((6,6) -> 1| 6 |<| 6 |(^)(6)(^ 2 )(6)) (h,t 2 ) 



dw x {i x )dw 2 {i 



2) 



&KI&I ieil Q - 1/2 l6l Q - 1/2 [i(£i+6)]p6] 



(1.6) 



and 



>l-(t) := 27rc a y dt 2 j 2 ^- 1 ((6,6)^l| 6 |<| 6 |(^)(6)(^2)(6))(ti,t 2 ) 
aWi(6)aW 2 (6) e 1 *^ 1 ^ 2 ) 



(1.7) 



Then 



As,t) = ^{( A+ ( t )- A+ ( s )) ~ (A-(t)-A-(s)) + (.4+( s ,i)-^M))} , 



where 



A+(s,t)-Ag(s,t) = l- J 



(1.8) 



+ 



Ci|<|&| 

( e H& _ e is6)gi^i | dWi(ii)dW 2 {&) 



(1.9) 

Two lines of computations show immediately that the variance of A% is 
always finite (essentially, considering e.g. Aq , because J^ 00 ^i+i J^°° ^f+i < 
oo, while the artificial infrared divergence at 6 = disappears when Tay- 
lor expanding e 1 ^ 1 — e 1 ^ 1 ). On the other hand, letting £ := 6 + 6 and 
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introducing an ultra-violet cut-off at = A 3> 1, one may see for in- 
stance A + (t) as an inverse random Fourier transform of the integral £ i— > 

dM g2 (S2) |^-g 2 | a -i/2| g2 | a -i/2 » whose variance diverges like / A ||| = 
0(A 1_4q ) or O(lnA) in the limit A — > oo as soon as a < 1/4. Note that the 
ultraviolet divergence is in the region \^\ ^> Later on, we shall in- 
troduce the field a, and this region will correspond to a splitting of the vertex 
J d<j> 1 {x)<j> 2 {x)a{x)dx = /d£id&d£<J(£i+6+£ = O)(J-00i)(&)^2(6)JV(O 
such that the momentum of the cr-field, £, is much lower than that of the two 
0-fields. Note also that, had we not operated the Fourier normal ordering, 
the above integral would have been infrared divergent near £2 = 0. 

It is apparent that the central role in this decomposition is played by the 
Fourier projection operator Z^li^^i^i) = J 7-1 (li^i^^i • T( . )) . Since 
A^j are obtained by Fourier projecting (-Bi(i) — Bi(s))4>2(s), or (^(t) — 
B2(s))4>i(t), which are perfectly well-defined products of continuous fields 
0, it was clear from the onset that these would be regular terms. Hence 
singularities come only from the one-time quantity A . (i), which does not 
split into a product of first-order integrals. 

In order to go one step further towards our final formulation, let us 
consider instead of the brute- force projection J^li^i^^i) the following pro- 
jection operators, which are more illuminating, and also analytically more 
robust, allowing for the extension of our theory to more general multiscale 
Gaussian fields as defined in section 2. 



Definition 1.3 (Fourier projections) Let tpi = ipi(x), i = 1,2 be Gaus- 
sian fields, and ipi = Yl'jLo their scale decomposition along the Fourier 
partition of unity (x J ) as * n Definition \2.10l Then one lets 

_. +00 

D+i^®^)^ Yl W 8^+2^'® ^2 (1-10) 

0<j<k<oo j=0 

and similarly 

j +00 

0<j<fc<oo j=0 

Clearly (D + + D~) (tpi ® tp 2 ) = ipi <8> ip2- 

8 apart from the spurious infra-red divergence (see above) 
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Thus one abandons - to the great benefit of notations - the above com- 
plicated expression for A^(t) in favour of the analytically equivalent 



A + (t) := D- 



-i +00 



0<j<fc<oo 



j=0 



A~(t) :=D-( / d<h(h)Mh)dh 



(1.12 

-1 +00 

£ E 9 « 



(ti)dti, 



0<j<fc<oo 



J=0 



(i x )dt 



1.2 Definition of the interaction 



i- 



(1.13) 



We recall that J dBi(t\) J* 1 ^.62^2) represents the area between the straight 
line connecting {B\{s), ^(s)) to (Bi(t), B2(t)) and the graph. If the curve 
turns right, resp. left, then the Levy area increases, resp. decreases. We 
have seen that A^ represents in some sense the singular part of the Levy 
area. So A + , resp. A" plays the role of a right-, resp. left-turning (singular) 
field. 

It is conceivable that B\,B% or 4>\,4>2 represent the idealized, strongly 
self-correlated motion in R 2 of a particle, which - although rotation-invariant 
- may not (probably as a consequence of a mechanical or electromagnetic 
rigidity due to the macroscopic dimension of the particle, or any other sim- 
ilar phenomenon) turn absolutely freely. A natural quantum field theoretic 
description of this rigidity phenomenon is to add an interaction Lagrangian 
of the form Ci n t = (dA^ 1 ) 2 . The fundamental intuition here is that the field 
B is in some sense a mesoscopic field, while A^ depends on microscopic 
details of the theory. This is explained in great accuracy in [21] , in a mathe- 
matical language. A. Lejay explains how the paths of B may be modified by 
inserting microscopic bubbles along the paths of B, resulting in the limit in 
paths which are indistinguishable from the original ones, while the Levy area 
has been corrected by an arbitrary amount. Hence one must search for an 
interaction which cures the ultra-violet divergences of the microscopic scale, 
without modifying the theory at mesoscopic scale. Understanding the meso- 
scopic scale as low-frequency (which is not really appropriate, since there 
are two reference scales here, instead of one), a natural candidate would 
be an interaction theory which is asymptotically free at large (mesoscopic 
scale) distances. The best-known example of that is probably the infra-red 
</> 4 -theory in 4 dimensions [9]; in that case, the coupling constant increases 
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indefinitely at short distances, and forces an ultra-violet cut-off. In our case, 
the coupling constant A will not flow, and the theory will be well-defined at 
all scales, which suggests a just renormalizable theory (or, in other terms, 
an integrated interaction which is homogeneous of degree 0). Since (dA^) 2 
is homogeneous of degree (4a — 2) in time, one shall use in fact a non-local 
interaction lagrangian, \c' a J J \ti - t 2 \~ 4a Ci n t((t>i, <h)(h, h), where 

Cua(<h,<h)(ti,ti) = X 2 {dA + (h)dA + (t 2 ) + dA-(h)dA-(t 2 )} , (1.14) 

which is positive for a < 1/4 since the kernel |ti —t2\~ 4a is locally integrable 
and positive definite. Thus the Gaussian measure is penalized by the sin- 

gular exponential weight ^ Cint ^'^ tl ' t2 ^ tl ~ t2 \~ 4adtldt2 . Equivalently, 
using the so-called Hubbard-Stratonovich transformation we introduce 
two independent exchange particle fields a± = o~±(t) with covariance ker- 
nel Ea±(s)a±(t) = c' a \s — t\~ 4a and rewrite (letting dn{4>), resp. dfi(a) be 
the Gaussian measure associated to <fi, resp. a = (cx + ,cr_)) the partition 
function Z := Z(X), 

Z (1.15) 

as 

Z := J e-f*^fa>tort®*dii(<l>)dvi(a), (1.16) 

where 

&nt(<l>i,<h,<r){t) = iX(dA + (t)a + (t)-dA-(t)a^(t)). (1.17) 



All of this is ill-defined mathematically since (1) a is a distribution- valued 
process and dA^ is not defined at all when a < 1/4; (2) one integrates over M. 
a translation- invariant quantity (note that <p\ , 02 , o~ are all stationary fields) . 



1.3 Towards a heuristic expression for the Levy area 

Assume the coupling parameter A is small enough. Perturbative quantum 
field theory suggests then to expand formally the exponential of the La- 
grangian and compute polynomial moments, tpi(xi) . . . t/j n (x n )e~ f £i™t(4>i,<p2,<?)(t)dt 
also called n-point functions and denoted by (ipi(xi) . . . ijj n (x n )) \, ipi = 



9 which is nothing else but an infinite-dimensional extension of the Fourier transform 
E[e iAX ] = e -" 2 ^ /2 for a random variable X ~ jV(0, a 2 ) 
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01,02,0+ or cr_, as \ J2n>o [4>i( x i) ■ ■■ r 4>n{x n ) (J C int ( ■ ;t)dt) n ]. Us- 

ing Wick's formula (recalled in §6.1), one may represent (ipi{x{) . . . ip n (x n )}\ 
as a sum over Feynman diagrams, ^r^(-O) where T ranges over a set of 
diagrams with n external legs, and A(T) £ C is the evaluation of the corre- 

§onding diagram (see section 4) . The Gaussian integration by parts formula 
yields a so-called Schwinger-Dyson identity, 



(dA ± (x)dA ± (y)) x 



1 



A 2 Z(A) 
1 



E 



-E 



Sa + (y) 5a+(x) 



- / Cint (01 ,<j>2,<T+)(t)dt 



■ f Cint (01 ,02 ,<T+ ){t)dt 



A 2 Z(A)~ l y ~ a +~ r ""'8a + (xy 
~ [-C£(x,y) + ((C- + V + )(x)(C- + V+)(y)) A 



(1.18) 



with Fourier transform 

(\T(dA±m\ 2 ) x = ^- ia [i - \^- ia (\(^ + m\ 2 )x\ . (i.i9) 




5-5, 



Figure 1: Bubble diagram. Bold lines are 0- fields, plain lines are cr-fields. 





Figure 2: First two terms of the bubble series. 

Consider (using a brute-force ultraviolet cut-off at momentum A) the 
term of lowest degree in A in the term between square brackets: it is essen- 
tially equal, up to a sign, to the evaluation of the half-amputated bubble 

10 an infinite-dimensional extension of the well-known formula for Gaussian vectors, 
E [d Xi F(X u . . .,X n )] = J2j C" _1 (i,i)E [XjF(Xi, . . .,X„)] if C is the covariance matrix of 
(Xi, . . . , X n ). 
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diagram, see Fig. dj 

_|£|l-4* . ( _ iA) 2 f A d ^{(E[\Ta + (0\ 2 }) 2 Efl^O^)! 2 ] E[|^ 2 (^-ei)| 2 ]} 

J\ti\<\t-ti\ L J 

= A 2 |e| 4Q - 1 / ^ikil^le-ar^-A^oo^A^A/lel) 1 - 40 , (1.20) 

J|£il<l£-£i 



hence a diverging positive quantity. (Using the Fourier truncation of Defi- 
nition [L3] only changes the constant K.) However, resumming formally the 
bubble series as in Fig. [2] yields - taking into account the possible insertion 
of cr_ -propagators between <r+-propagators - 



J_|£|l-4a 



1 



1 + K>\*(AM) 



l-Aa 



i x _ 4q ^(A/iei) 1 - 4 " 



A 2ISI l + K'A^A/lel) 1 - 4 " 

-+A->oc ^ICI 1 - 4 "- (1-21) 



Thus the bare u-propagator | ^ | has been replaced with the renormal- 
ized propagator ic-\-KX^A^~ ia ' which vanishes in the limit A — y oo. In 
physical terms, the interaction in |g|i-4a has been totally screened by an 
infinite mass term K' A 2 A 1_4q . More complicated diagrams contributing to 
(|(J r <7+)(£)| 2 ),\, and involving internal cr-links, also vanish when A — >• oo. 
There remains simply: 

(l^ ± (0l 2 l>A = ^|€|- 1 - 4tt . (1-22) 
As for the mixed term (dA ± (x)dA T (y)}\, its Fourier transform is given 
by ^|£| 1-4a — i + K"\' 2 (A/\£\) i - ia ' wnere K" < K' due to the constraints 
on the scales for bubbles of mixed type with one a+- and one <7_ -leg 0, 
which vanishes in the limit A — > oo (note the disappearance of the factor 1 
compared to eq. (11 .21 f) . due to the fact that Ecr + (x)<7~ (y) = 0). Thus the co- 
variance of the two-component cr-field has been renormalized to r^rr^kj^f^ , 

where b p is a two-by-two positive matrix with eigenvalues ~ A 2 M p ( 1_4Q ) . 
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Namely, coupling (fa, fa) simultaneously to a+ and a_ leaves only the scale-diagonal 



projected field D + D~ (fa ® fa) = \ I] fa ^ ^ 
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Using eq. (jl.8p . one obtains: 

(27rc a ) 2 (A(s,t) 2 ) x = (\A + (t)-A + (s)\ 2 ) x + (\A-(t)-A-(s)\ 2 } x 

+E\A+(s,t)-Ag(s,t)\ 2 

= 1 j\l-co S (t-s)0\t\- 1 - 4a dt + E\A+(s,t)-AQ(s, 
= (±K 1+ K 2 )\t-s\ ia (1 
for some constants K\ , if 2 ■ 




Figure 3: Higher connected moments of the Levy area. 

A simple power-counting argument (see subsection 4.1) shows that the 
overall degree of divergence of a connected diagram with 2n external <7-legs 
is 1 — Ana. For n > 2, this is < 1 — 8a < since a > g by hypothesis, so 
such diagrams are convergent. By the above arguments, there remain only 
the connected diagrams in the limit A — > 00, see Fig. [3l whose evaluation 
is independent of A. Then these may be shown by standard multi-scale 
arguments for single Feynman diagrams (see e.g. [39]) to be of order K n , 
which implies that the generating series for connected moments has a non- 
zero radius of convergence. 

On the other hand, the law of the field <fi is left unchanged by the inter- 
action. Namely, all non-trivial diagrams contributing e.g. to (4>i(x)(f)2(x)) \ 
involve internal u-links which (as previously "shown") vanish in the limit 
A — > 00. 

On the whole, this is the content of Theorem 0.1. 

The art of constructive field theory is to make the previous speculations 
rigorous. 
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2 Multiscale Gaussian fields 



2.1 Scale decompositions 

We fix some constant M > 1. The next definition is borrowed from [31] . 

Definition 2.1 (Fourier partition of unity) Let x° '■ K — > [0,1] be an 
even, C°° function such that X°\^_i ji = 1 an ^ SU PPX° C [— M, M], and 

X^O^X^M-^-X^M 1 "^), i>l (2-1) 

so that suppx j C [M 3 '"~ 1 ,M J ' +1 ] U [-M^ +1 , — M J ' -1 ]. 

TTie (Xj)j>o define a C°° partition of unity, namely, Xlj=o Xj = 1- 

Note that suppxj H suppxj' has empty interior if |j — j'| = 2, and is 
empty if \ j - j'\ > 3. 

Definition 2.2 (M-adic intervals) Let W := {[kM~i , (k + l)M~ j ),k E 
j > be the set of M-adic intervals of scale j, and D := ttljxjB^ 
i/ie disjoint sum of these sets over all scales. The set D is a tree with 
links ( called: inclusion links ) connecting each interval A E W to the unique 
interval A' E W~ l such that A C A' (see Definition \3. 7| below, and left part 
of Fig. 2] subsection 3.3). An element of W is usually denoted by Ai , 
or simply A if no confusion may arise. The volume \A 3 \ is simply M~ 3 . If 
A E D- 7 , f/ien one denotes by j(A) = j £/te sca/e of A. 

If x GM, then x belongs to a single M-adic interval of scale j, denoted 
by A j x . 

If A J E B J ; i/jen the set of intervals A E l±l/ l< jB' 1 suc/i i/iai A lies below 
A j , i.e. AD A j , is denoted by (A- 7 ')^. 

We denote by d J (A, A'), A, A' e D J , £/te distance in terms of number of 
M-adic intervals of scale j between A and A', namely, 

d j ([kM~ j , (k + 1)M- J ), [k'M~ j , (k' + l)M~ j )) = \k' - k\. (2.2) 

By extension, one may also define the d? -distance of two points or two 
M-adic intervals of scale j' > j, namely, d^(x,y) = M^\x — y\ and 

d j (A, A') := M j - j 'd j '(A, A'), A, A' E B J ". (2.3) 

Remark. It is preferable not to define dP(A,A') for A, A' E B J ' with 
j' < j, since d J (x,x'), x E A,x' E A' depends strongly on the choice of the 
points x, x' then. 

The definition extends in a natural way to a D-dimensional setting by 
decomposing M. D into a disjoint union of hyper cubes of size side M~K 
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Definition 2.3 (multiscale decomposition) (i) Lettp = tp{x) be some 
integrable function or distribution in S' . Then the multiscale decom- 
position of tp is tp = Y2j>o i* 3 > w here 

P(x) := T- 1 {$ ^ xHO^W) (x) = ^=J ^(0(^)(0 di. 

(2.4) 

(ii) (low-momentum field) Let tp~* k := Y^]j=o' l l )J ■ 
(Hi) (high-momentum field) Let tp k ^ := ^2^P^.tp- 

Now comes a general remark. Let / = J-f(£) be some function with 
support in |f | < M' such that \T{f '){£,)] = is bounded. Then 

\F- 1 f(x)-F- l f{y)\<\x-y\ / \C\\f(0\ d^ < K ■ M^\x - y\, (2.5) 

Jo 

so ip 3 (x) or ip~*3(x) varies slowly inside intervals of scale k if k > j. 
Hence it makes sense in first approximation to consider tp 3 {x) or tp~* 3 {x) to 
be approximately equal to the averaged, locally constant function 

iPl(x):= £ K&arhr f ^(y)dy = [ ^(y)dy, (2.6) 
or similarly for the low-momentum field tp^ 1 

1>?v S {*)-= E ^ TX7T / ^ j (y)dy = -±j- [ .^ j (y)dy. (2.7) 

Summing the ^ u over j would give a new function ^j>o ^wu ( x ) which is 
a sort of "naive" wavelet expansion of the original function tp; with a little 
extra care, one could arrange that the two functions be equal, but we shall 
not need to do so. 

Conversely, if ip is a 'reasonable' random Gaussian field, then the covari- 
ance {tp J (x)tp J (y)) — or, more generally, {tp 3 ^ (xjtp^ (y)) - is usually small 
if the corresponding M-adic intervals are far apart, i.e. if d 3 (A J x , Ay) S> 1. 

These two remarks may be made precise in the case when tp is a multi- 
scale Gaussian field, see Definition 12.41 below. Then, for any scale j, 
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• the field ip^ J may be decomposed into the sum of the locally averaged 
field at scale j, namely, ipau{x), and a secondary field, denoted by 
S^i/j^ 1 , whose low momentum components of scale h < j decrease like 
Af-TU'-' 1 ) for some 7 > 0, see Lemma 12.61 and Corollary 12. 71 For 
reasons explained below, it is customary to use the nickname of spring 
factor for a decrease factor of the type M~^~ h \ 7 > 0. 

• since the covariance decreases with the distance in terms of number 
of M-adic intervals of scale j, it makes sense to try and find some 
expansion of the functional C(tp) in which the field values over far 
enough M-adic intervals have been made independent. This is called 
a cluster expansion (see section 2). 

Definition 2.4 (multiscale Gaussian field) A multiscale Gaussian field 
with scaling dimension [3 < 1 is a field if) = ip(x) such that, for every r > 
and r, t' = 0, 1, 2, . . 

\(d^( X )dr $>(y))\ < K T y, r -—— (2.8) 

(1 + M3\x — y\) T 
with some constant K TT i T depending only on t,t' and r. 

Remark. A multiscale Gaussian field ip with scaling dimension (3 £ 
(0,1) has almost surely /3~-H61der paths. Namely, K\ipi(x) — ip 3 ' (y)\ 2 = 

f% dzdz'(dipi(z)di/ii(z r )) is bounded by K(x - y) 2 M^- 2 ^ if \x - y\ < 
M~ j , and by 

r\x-v\ r\x-y\ M (i-2f3)j 

2j dz'\(iP(0)diP(z'))\ <Kj o dz' {1 + Mjz/)2 < ICM-W (2.9) 

otherwise. Summing over j yields 



E\iP(x) - ^(y)\ 2 < K \{x - yf M ^~ 2 ^ j + M ~ Wj 

\ J<— log |cc— y| i>-log|cc-y| 

< K'\x-y\ 2 P, (2.10) 

and one concludes by using the Kolmogorov-Centsov lemma. 

As we shall see in the next paragraph, fractional Brownian motion with 
Hurst index a is the paramount example of multiscale Gaussian field with 
scaling dimension [3 = a £ (0, 1). 
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Definition 2.5 (averaged and secondary fields) Choose some scale k > 
0. Let f : R — > R 6e a function, typically, f(x) = ip~^ k {x) oripi(x) for some 
j < k. Then: 

(i) the averaged field at scale k is the locally constant field 

/ f(y)dy. (2.11) 

It is more convenient to consider f av as some function on D k , still 
denoted by f , 

f(A k ):=j^ r J^f(y)dy, A k €B k , (2.12) 

so that f av {x) = f{A^,). This notation fixes unambiguously the scale 
of the average. 

(ii) the secondary field at scale k is the difference between the original field 
and the averaged field at scale k, namely, 

6 k f(x):=f(x)-f(A k ). (2.13) 

These definitions imply the following easy Lemma: 
Lemma 2.6 

S k f(x) = J^rJ Ak /'(«)*>) du = dvf'(v)5 k (x; v) (2.14) 
where 

5 k (x;v) := {(u-infA^l^ + ^-supA*)!^} G [-1,1] (2.15) 

is a signed distance from v to the boundary of the interval A k = [inf A^, sup A^) 
measured in terms of the rescaled d k -distance. 

Proof. Straightforward. □ 

Corollary 2.7 Let k,k' > j. Assume ip = ip(x) is a multiscale Gaussian 
field with scaling dimension f3 < 1. Then 

\{6 k ^ j (x)5 k '^ j (y))\ < K r - — M-^M-^' -^(2.16) 

i\ v v ) v \yj/\ - r (i + M3\x-y\) r v ' 

= K r M~^ k+k r- ^ . (2.17) 

(l + Mi\x-y\) r V ' 
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Proof. Straightforward. □ 
Eq. (|2.16p . resp. (|2.17p emphasizes the "spring-factor" M~( k ~i\ resp. 
M~^~~^ k ~i> gained for each secondary field with respect to the covariance 
of a multiscale Gaussian field with scaling dimension /3 at scale j, resp. 
at scale k. Had one considered directly \(ifr> (x)^ (y))\, the spring factor - 
called rescaling spring factor - in (|2.17p would have been simply M^ k ~^ , 
see introduction to §6.1.2. 
Remarks. 

1. The same spring factors appear in D dimensions. It is sometimes 
useful to take a cleverer definition of the secondary field by replacing 
the simple average /(A^) with a wavelet component of /, where the 
wavelets have vanishing first moments up to order r > 1. This allows 
one to enhance the spring-factor from M~^ k ~^ to M~( T+l ^ k ~i\ resp. 
from M-^-flfc-fi to M- ( - T+1 ~^ k - j \ 

2. The separation of low- momentum fields into a sum (field average) + (secondary 
field) is required only for fields with scaling dimension /3 > —D/2, and 

is not performed otherwise (see explanation after Definition 13.91 and in 
subsection 6.1). 

Consider conversely a high-momentum field ipi(x) : j > h with x G H> h . 
Typically (see section 3), ipi(x) = ^(x^h), ^{x' lh ) or ^ ( x A h ,t) i x e A ft 
coming from the scale h horizontal cluster expansion or a scale h i-derivation 
in an interval A h G H) h , and one must integrate over x G A h . Then ipi(x') 
and ipi(x") are almost decorrelated if x',x" G A h but d^(x,x') 3> 1. Hence 
it makes sense to restrict ipi over each sub-interval A- 7 C A h of scale j: 

Definition 2.8 (restriction of high-momentum fields) Let A h G H> h 

and j > h. Then the high-momentum field ipi{x), x G A h , splits into 

ifj(x)= Yl Res^Tp j (x), x€B h (2.18) 

AJ'Gl!«',AJCA' 1 

where Res^jifr' (x) := i x eAo^ ( x )- 

2.2 Multiscale Gaussian fields in one dimension 

We introduce here more specifically the infra-red divergent stationary field 
<j> associated to fBm B, and the fields a = a± conjugate to the turning fields 
A^. In this paragraph D = 1. All fields come implicitly with an ultra-violet 
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cut-off at scale p, so that ip = (ft, a should be understood as ip~ >p , Tp J ^ as 
ipi^P, and so on. 

Definition 2.9 (Harmonizable representation of fBm) LetW "(£)>£ ^ 
M lie a complex Brownian motion such that W(— £) = W(£), and 

i r+°° e 1 ^ - 1 i 
B t := (2Trc a )-2 / —^\(:\2—dW(0, teR. (2.19) 

J-oo K 

The field i 6 K is called fractional Brownian motion\^ . Its paths are 
almost surely a - Holder, i.e. (a— e)-H61der for every e > 0. It has dependent 
but identically distributed (or in other words, stationary) increments Bf — B s . 
In order to gain translation invariance, we shall rather use the closely related 
stationary process 

r+oo it£ 

/ -fI^"^?), (2.20) 

- with covariance 

<0(a#(y)) = | ^(-»)^l_de (2.21) 

- which is infrared divergent. However, as already discussed in section 1, 
the increments <p(t) — (j>(s) are well-defined for any pair of variables (s,t). 

Definition 2.10 1. Let 

Cjfoy) := J e ^)|ig-^, j E N. (2.22) 

Then := X^'eN ^6 * s ^ e covariance of the field <f>. We denote by 
4> '■= X^'eN 4? the corresponding multiscale decomposition of the field 
(j) into independent components ^, j 6 N. 

2. Let cjf^i = 5^1=0^ anc ^ = <fc'^ p = J2k=j^ k - The covariance 
of 4>, resp. (jy* 3 , resp. (j> > ~* is D(x')C ( j ) , resp. D(x ->J ")(7^, resp. 
D( X j -+)C4, where X ~* j := E{= X h , X j ~> = ELjX k - 

Remark. Note that this multiscale decomposition does not coincide 
exactly with that of Definition 12.31 With this slightly modified definition, 
the (fP have been made independent. 

12 The constant c a is conventionally chosen so that ~E(B t — B s ) 2 = \t — s\' 2a . 
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Lemma 2.11 The stationary field (j) associated to fBm is a multiscale Gaus- 
sian field with scaling dimension a. 



Proof. A simple scaling of the variable of integration yields 

C$(x,y) = M- 2 ^- 1 ) J e^-^-y^'j^xHOde (2-23) 

withsuppx 1 C [1,M 2 ]U[-M 2 ,-1] bounded away from 0, hence \C^(x, y)\ < 
M~ 2a i , and more precisely (by integrating by parts n times) \C^(x,y)\ = 
O (M- 2a i(Mj\x - y\)~ n ) when M j \x-y\ -> oo. The bound for {d T & (x)d T ' & (y)) 
is obtained in the same way (simply multiply by |£| T+r in Fourier coordi- 
nates) . □ 

We may now define the field a. As indicated in the Introduction and 
in section 1, the Fourier transform of its bare covariance (<7j<Tj/), = ± 

is ■ O n the other hand, the renormalized covariance (up to the 

overlap between the supports of the Fourier multipliers x 3 ) 1S essentially 

mi d„Tj , v^=Fg .,■ Dm i where the scale j mass counterterm W for dia- 

l^i *"id+X, 3 - = o^X J (?) 

grams of scale j with two low-momentum external legs cr - *^ -1 ) - a two- 
by-two, positive matrix - is defined inductively (see subsection 3.5, section 
4 and subsection 6.3) and shown to be of order A 2 M^ l ~ Aa \ For technical 
reasons one chooses to retain in the covariance of a only a simplified version 
of the counterterm (essentially the term of highest scale p), which is of the 
same order as the sum of all mass counterterms. 

Definition 2.12 (i) Let a be the stationary two-component massive Gaus- 
sian field with covariance kernel 

r+oa p i£(x-j/) 

c °' (x - y)i =L if (2 - 24) 

(ii) Decompose C^ p into ^2j >Q Ci and a into a sum of independent fields, 
X^'>o cr ' J as in Definition \2.10\ by setting 



CUx,y) = f e^t'-w) df. (2.25) 
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Lemma 2.13 a is a multiscale Gaussian field with scaling dimension —2a. 
More precisely, 



. A/f (r+r'-4a)j lirj(l-4a) 

(d T o-i(x)d T a\y)) <K TT , r - — - ^ . inf(l, ). (2.26) 

\ v i vy;/ - (1 -)- A13 \x — y\) r bP ' ' 

For p large enough, inf(l, MJ( ^~ 4a) ) w \-2 M ~( P -j)(i-^) _ Thug ff van _ 
ishes in the limit p — > oo because of the infinite mass counterterm. 

Proof. Same as for Lemma 12.111 Note that the rescaled denominator 
|g|i- to+b pM (g L~i)(i-4 a) is bounded by inf(l, MJ( ^ 4a) ). □ 

Note the following to conclude: 

- multiscale Gaussian fields with positive scaling dimension (as (j) for 
instance) are ultra-violet convergent and define continuous fields, but may be 
infra-red divergent, necessitating some infra-red cut-off (however, if /3 < 1, 
the increments 4>t — 4>s are well-defined); 

- multiscale Gaussian fields with negative scaling dimension (as in the 
case of the field a, or of most fields in quantum field theory) are ultra-violet 
divergent and define distribution-valued fields whose covariance kernel is 
singular on the diagonal. On the other hand, they are infra-red convergent. 



3 Cluster expansions: an outline 

The general aim of horizontal and vertical cluster expansions has already 
been explained in the Introduction and in §2.1. The (non restricted) hor- 
izontal cluster expansion has been given by D. Brydges and T. Kennedy a 
beautiful combinatorial structure in terms of forests of intervals (see §3.1 
and 3.2). The vertical or momentum- decoupling expansion, in terms of t- 
parameters, is somewhat looser, relying on a Taylor expansion to some order 
ta in each interval A. Putting together these two expansions, one obtains 
so-called polymers (see §3.3). Polymers with too few external legs must still 
be renormalized by resumming the local part of diverging graphs into an ex- 
ponential making up an interaction counterterm (see section 4 for detailed 
explanations); a Mayer expansion makes it possible to get rid of their non- 
overlap constraints with the other polymers, and finally to divide out the 
polymers with no external legs (called: vacuum polymers) when computing 
n-point functions (see §3.4). All these provide series depending on the small 
parameter A which will be shown to converge in section 6. 
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Summarizing, the general expansion scheme which one must have in 
mind is made up of a sequence of transformations for each scale, starting 
from highest scale p, namely 

Horizontal cluster expansion at scale p > Vertical cluster expansion at scale p ^ separation of local 

part of diverging graphs of lowest scale p > Mayer expansion at scale p ^ resummation of local parts of 

diverging graphs of lowest scale p, also called : rcnormalization stage of scale p ^ > . . . > ^Horizontal cluster 

expansion at scale j > Vertical cluster expansion at scale j > separation of local part of diverging graphs of 

lowest scale j > Mayer expansion at scale j > resummation of local parts of diverging graphs of lowest scale 

*)-... 

down to lowest scale j = 0. 

3.1 The general Brydges-Kennedy formula 

Let us start with the following general definition. 

Definition 3.1 Let: 

(i) O be an arbitrary set, whose elements are called objects; 

(ii) L{0) be the set of links of the total graph associated to O, or in other 
words, the set of pairs of objects, so that I G L{(D) is represented as a 
pair, i ~ {o£, o' £ } C O, og / o' e ; 

(Hi) [0,1] L ^ := {z = (^)£gL(O))0 < z# < 1} be the convex set of link 
weakenings of O; 

(iv) F(0) be the set of forests connecting (some, not necessarily all) ver- 
tices of O. A typical element of F{0) is denoted by F, and its set of 
links by L(F) C L(O). 

Assume finally that link weakenings have been made to act in some 
smooth way on Z (a functional depending on some extra parameters), thus 
defining a C°° link-weakened functional Z : L(O) — > M, z = (ze)£ e L(o) ^ 
Z(z) on the set of pairs of objects, still denoted by Z by a slight abuse of nota- 
tion, such that Z (the original functional) is equal to Z(l, . . . , 1) = Z{1 L ^). 

Definition 3.2 (one step of BK expansion) The Brydges-Kennedy de- 
coupling expansion consists in the following steps: 
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(i) Taylor- expand Z(l, . . . , 1) with respect to the parameters (zi) simulta- 
neously, namely, 

Z(l,...,l) = Z((z e ) ieL{O )=0) (3-1) 

+ Yl / dw ^- 9 H 1 Z (( z i)eeL(0) = wi) ■ 
e 1 eL(O) j0 

(ii) Choose some link i\ in the above sum. Draw a link of strength w\ 
between oe 1 and o'g and consider the new set 0\ made up of the simple 
objects o G O \ {o^, o' ti } and of the composite object {(o^, o' ti )}, with 
set of links L(Oi) = L(0) \ Then d Zt ^Z (z h = w\, (^)^l(Oi)) 

must now be considered as a functional of (ze)eeL(o 1 )- 

When iterating this procedure, composite objects grow up to be trees. 
This leads to the following result: 

Proposition 3.3 (Brydges-Kennedy or BK1 cluster formula) Z(l, . . . , 1) 
may be computed as an integral over weakening parameters W£ € [0, 1], where 
£ does not range over the links of the total graph, but, more restrictively, over 
the links of a forest on O : 



z(i,...,i)= £ ( n fdwA (( n M 

¥eF(o) \ieL(¥) Jo j \ \eeL{¥) e J j 

where zg{w\l G L(O) is the infimum of the for I' running over the 
unique path from 0£ to o'^ if oi and o'^ are connected by ¥, and z^{w) = 
otherwise. 

If it is not desirable to make a cluster expansion with respect to the 
links between certain objects (of type 2 in Proposition 13.41 below) . then it is 
sufficient to consider all these objects as belonging to the same composite 
object. This yields: 

Proposition 3.4 (restricted 2-type cluster or BK2 formula) Assume 
O = 0\ II 02- Choose as initial object an object o\ G 0\ of type 1, and stop 
the Brydges-Kennedy expansion as soon as a link to an object of type 2 has 
appeared. Then choose a new object of type 1, and so on. This leads to a 
restricted expansion, for which only the link variables Z£, with £ Oi x Oi, 
have been weakened. The following closed formula holds. Let J- res {0) be the 
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set of forests ¥ on O, each component of which is (i) either a tree of objects 
of type 1, called unrooted tree; (ii)or a rooted tree such that only the root 
is of type 2. Then 

zd,. ..,D= £ n n l 

(3-3) 

where Sg(w) is either or the minimum of the w-variables running along 
the unique path in ¥ from oe to o'^, and ¥ is the forest obtained from ¥ by 
merging all roots of¥ into a single vertex. 

This restricted cluster expansion will be useful for the Mayer expansion 
(see section 3.4). 



3.2 Single scale cluster expansion 

Definition 3.5 (i) A horizontal cluster forest of level p' < p, associated 
to a d-dimensional vector Gaussian field ip = (ipi(x), . . . ,ipd(x)), is a 
finite number of M-adic intervals in © p , seen as vertices, connected 
by links, without loops. Any link I connects Ag to (A^, £ O p ), 
bears a double index, (ie,i'g) £ {1, ...,d} x {l,...,d}, and may be 
represented as two half-propagators or simply half-segments put end to 
end, one starting from Ag with index ig, and the other starting from 
A^ with index i'«. 

The set of all horizontal cluster forests of level p' is denoted by J- p . If 
¥ p ' G F p ' , then L(¥ p ') is its set of links. 

(ii) A horizontal cluster tree is a connected horizontal cluster forest. Any 
horizontal cluster forest decomposes into a product of cluster trees 
which are its connected components. 

(Hi) If there exists a link between A and A' (A, A' G D p ) then we shall 
write A ~ FP / A', or simply (if no ambiguity may arise) A ~ A'. 

The following result is an easy consequence of Proposition 13.31 see pQ. 

Proposition 3.6 (single-scale cluster expansion) Let 

Z v {\):= [ e~fv WW** dfi y,-) t (3.4) 
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where dp is the Gaussian measure associated to a Gaussian field with d > 
1 components, ip{x) = (i/ji(x), . . . , ip^x)), with covariance matrix C = 
C(i,x;j,y), and C is a local functional. 
Choose p' < p. 

Let dp, s (ip), s = (•SA,A')(A,A')eBp'xDp' e I ' 1 ] such that S A,A' = sa'.a 
and sa,a = 1; be the Gaussian measure with covariance kernel (if positive 
definite ) 

C s (i, x; i',x') = s Ax ,A x ' C (h x; i', x') (3.5) 

if A x 3 x, resp. A x > B x' (A x ,A x i G D p ') are the intervals of size M~ p ' 
containing x, resp. y. 
Then 



Z V (\):= Y, 
fp' eFp' 



\\ dw e dx £ / dx' e C s(w) (ii,x e ;i' e ,x' e ) 

_rZ , s Jo JA e J A' 



r///. s(u „(r)TIo^(e-^ £ W^ da: ) (3.6) 

where: 



(3.7) 

s(w) = (sa,A'( w ))a A'ei/' S A,A'{ W )> A / A' being the infimum of the W£ 
for I running over the unique path from A to A' in ¥ p ' if A ~ FP / A', and 
S A,A' i w ) = else. 

Note that C s ( w j is positive-definite with this definition [I], as a convex 
sum of evidently positive-definite kernels. 



3.3 Multi-scale cluster expansion 

As explained in the Introduction, cluster expansion has two main objectives. 
The first one is to express the partition function as a sum over quantities 
depending essentially on a finite number of degrees of freedom. For a given 
scale j, the horizontal cluster expansion perfectly meets this aim. The second 
one is to get rid of the invariance by translation, which necessarily produces 
divergent quantities. 

Multi-scale cluster expansion fulfills this program by building induc- 
tively, starting from the highest scale p and going down to scale 0, connected 
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multi-scale clusters called polymers. These are finite, connected subsets of 
D, extended over several scales, and containing at least one fixed interval A-? 
at the bottom which breaks invariance by translation. Like the horizontal 
cluster expansion, they are obtained by a Taylor expansion with respect to 
some t-parameters (^a)a depending on the interval. 

Let us give a precise definition of such an object (see left part of Fig. 2]). 




P-4 | 

9 



— p-5 



Figure 4: Polymer. A horizontal bold line below an interval means that it 
is not connected to the intervals below it. 

Definition 3.7 (i) A polymer of scale j (typically denoted by ¥ J ) is a 
tree of M-adic intervals of scale j, i.e. a connected element of J- 3 . 
The set of all polymers of scale j is denoted by V J . 

(ii) A polymer down to scale j (typically denoted by W 3 ~*) is a connected 
graph connecting M-adic intervals of scale k = j, j + 1, . . . , p, whose 
links are of two types: 

- horizontal links connecting M-adic intervals of the same scale; the 
restriction ofWi~* to any given scale, P J_> n B fc , k > j, is required to 
be a disjoint union of polymers of scale k, or simply in other words, 
an element of J- k ; furthermore, P^nlP' is assumed to be non-empty; 

- vertical links or more explicitly inclusion links connecting an M-adic 
interval A E B fc , k > j to the unique interval A' E D fc_1 below A, 
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i.e. such that A C A'. These may be multiple links, with degree of 
multiplicity ta. 

A polymer F 3 ^ is also allowed to have an external structure, charac- 
terized by a subset A J C P J_> n W of intervals of scale j, such that 
each A 3 G A J is connected below by an external inclusion link to ED 3 ' -1 . 

f m ) The horizontal skeleton of a polynomial T 3 ^ is the disjoint union of the 
single-scale cluster forests P^nD*, k > j (all vertical links removed). 
One denotes by V 3 ^ the set of polymers down to scale j, and V J Next C 

the subset of polymers with N ext := ^A^ePnAi r A external links. 
In particular, polymers inV 3 ^ , without external links, are called vacuum 
polymers. 

The easiest example is that of so-called full-inclusion polymers F 3 ^ G 
V 3 ^ containing all inclusion links. Such polymers may be obtained from 
arbitrary multiscale horizontal cluster forests F 3 ^ = (F 3 , F 3+1 , . . . , F p ) by 
linking all pairs (A, A'), A, A' G F 3 ^ such that A' is the unique interval 
lying below A. 

In general, if A fc C A fe_1 , A k G V 3 ~^, then there is an inclusion link 
from A k to A k ~ 1 inside V 3 ^ (by definition) if and only if T^k / 0. 

The integer ta corresponds in the momentum-decoupling expansion to 
the number of derivatives d/dt&. 

Whereas horizontal links are built in by independent horizontal cluster 
expansions at each scale, the construction of vertical links depends on a 
procedure called momentum decoupling (or sometimes vertical cluster) ex- 
pansion, which we now set about to describe. 

Let us first give some definitions. 

Assume positive numbers t? x , < j < p, depending on x but locally con- 
stant in each M-adic interval A 3 G U> 3 , have been defined. Let us introduce 
operators T^ k (0 < k < p, x G R) acting on the low-momentum components 
of the fields at the point x, i.e. on <f> ' (a?) and cr 3 (x), j < k, by 

(T^) 3 (x) = T^ 3 (x) = t k x t k ~ x . . . 4+V (*), k > j, (3.8) 

(T~> j il>y(x) = T^ j ift 3 \x) = ^(x), (3.9) 

where ip = (ft, deft or a. 

The shorthand (T^ k ip) 3 emphasizes the idea that T~* k is not simply a 
multiplication by some product of t-variables, but an operator acting diag- 
onally on the whole field ift, or rather on ift~* k (seen as a vector). 
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In other words, writing t J x = t A j if x € A- 7 , so that i is a real- valued func- 
tion on W+P = tt) <j<pBJ, one has for k > j: (T^ k ipy(x) = Uk'=j+i ' 
ip 3 (x). 

We shall also use the following notation: 

(TV>r fe := ^(T^V) j - (3-10) 

This weakened field, depending on the reference scale k, will be called 
the dressed low-momentum field at scale k. Separating the ^-component 
and the i fc -variables from the others yields equivalently: 

(TV>r fc (x) = ^ h (x) + t k x ^(T^-^'Oe). (3.11) 

j<k 

Dressed interactions may be built quite generally out of dressed low- 
momentum fields. Let us give a general definition. 

Definition 3.8 (dressed interaction) (i) Let C'^' p {'4>){x) := X Kq tpi q (x) = 
A K? riiG/g ^i( x ) be some arbitrary local functional built out of a prod- 
uct of fields. Then the momentum decoupling of C q , or simply dressed 
interaction, is the following quantity, 

V(.;i)(s) := A** | J] WW + 5> " WiT^r^'^Xx). 

yi&Iq p'<p iElq 

(3.12) 

(ii) (generalization to scale- dependent case) 

Let Ct p mx) := (\T q E( ji)ieIq YlieijH*) be a local func- 

tional with coupling constant \ p and coefficients Kg depending on 
the scales. Then one defines 

V( . ;t)(x) := (A')"' Yl K t ] \{{T^r{x) 

0<(ji) ie i q <p i£l q 

+ ^ (A P'-i)^(l _ ( t P')l^l) Y, HiT^'-V^ix). 

p'<p o<(ji) ieIq < P '-i iei q 

(3.13) 

where X p ' , p' < p are the renormalized coupling constants. 
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Summing up in general the contribution of the different interaction terms, 
£int '■= Ylq=i Kq£q, this defines a dressed interaction C iri £( . ;t) and a 
dressed partition function 

Z^ p {\-t) = [ e -Jv^( • VW^dniip). (3.14) 



For the moment, vertical links have not been constructed, and nothing 
prevents a priori the different horizontal clusters to move freely in space one 
with respect to the other, nor on the contrary to generate via Wick pairings 
infinite multi-scale clusters. As we shall now see, the momentum-decoupling 
expansion provides the mechanism responsible both for the translation- 
invariance breaking and the separation of polymers. 

Since this is an inductive procedure, let us first consider the result of 
horizontal cluster expansion at highest scale p. It may be expressed as a 
sum of monomials split over each connected component T p , . . . ,T P P of ¥ p , 
with generic term (called G-monomial or simply monomial) 



>(p-l) 



(3.15) 



for some (possibly empty) index subsets Ie, Ji, J' e C {1, . . . , d}, multiplied 
by a product of propagators as in eq. (|3.6p ; the f-variables dressing low- 
momentum components have been written explicitly as in eq. (13. lip . 

Let us draw an oriented, downward dashed line from A p € T£ to some 
A p ' g (A^) u , p' < p below A p if G p c contains some low-momentum field 
tjj p (xi) with xe G A p . Either A p ' G" F p ', or A p ' belongs to some connected 
component T p , of F p ' . In the latter case, one has attached T p to some cluster 
tree below, Tj , by the inclusion constraint A p C A p ' , which prevents T p 
from moving freely in space with respect to Tj. It may also happen that 
G p c contains no low momentum field component ip p , p' < p, so that T p looks 
isolated; unfortunately, nothing prevents some horizontal cluster expansion 
at a lower scale p' from generating some high-momentum field component 
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ip p {x), x G A p ' G (A^) u , which may be represented as a reversed upward 
dashed line from A p to A p . 

In order to have an effective mechanism of separation of scales, we shall 
make a Taylor expansion to order 1 with respect to the (i A ) AgT p-variables 
of the product (G-monomial) x ((i A ) A6T p-dependent part of the dressed 
interaction), G p := G p e ^ C mt (x)dx ^ name jy (splitting G p into a product 
II a = 11a G p A e~ f& Cint ^ dx over the fields located in each interval A of 
scale p in T p ) 

G p A (t p A = 1) = G A ((i A ) AeT £ = 0) + jf 1 <^G A (t A ), (3.16) 

thus producing a new set of monomials multiplied by the interaction. 

Setting all (i A ) Ae Tg to zero has the effect, see eq. (|3.12p and (|3. 15[) . 
of killing in the interaction - and hence in G p which is a derivative of the 
interaction - all mixed terms containing both ip P (x) and ip p (x'), with x G 
AP C W, x' £ Ap' C Op', Ap' d Ap as above. Hence, in all the terms 
G p ((t A ) Ae jp = 0), Tc has been effectively isolated from all other intervals 
in ID; it constitutes a (single-scale) isolated polymer. Thus, by letting t& = 0, 
one cuts all dashed lines crossing from A (or above A in general) to A^, 
and sets up a wall between A and A^. 

On the contrary, derivating with respect to t A , A G T p , produces neces- 
sarily low-momentum components, and hence vertical links materialized by 
dashed lines as above. In this case, one links A p to A p ~ l by a full line, sig- 
nifying that all corresponding monomials will contain some low-momentum 
field; this implies in turn the existence of a downward dashed line or wire 
connecting some A p G T p to some A p G (A p )^ as above (note however that 
there isn't necessarily a dashed line from A p to A p_1 ). 

Due to the necessity of renormalization (see section 4) and to the domi- 
nation problem (see §6.2), one may need to Taylor expand to higher order, 
up to order N ext ^ m ax + 0(n(A)). One obtains thus a polymer with a certain 
number of external legs per interval. Choosing the Taylor integral rest term 
for some A leads to a polymer with > N ex t }ma x external legs, which does not 
need to be renormalized. On the other hand, a polymer whose total number 
of external legs is < N e xt,max requires renormalization (see section 4) . 

Let us emphasize at this point that high-momentum fields have two scales 
attached to them: j and k for a field ipf(x), x G A- 7 (k > j), produced by the 
horizontal/ vertical cluster expansion at scale j; but low-momentum fields tp^ 
have three scales. There are in fact two cases: 
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(i) If fa < — D/2 then tpi is not separated into a sum (low- momentum 
field average) + (secondary field). The genesis of (contrary to the 
high-momentum case) is actually a process which may not be under- 
stood apart from the multi-scale cluster expansion. At their produc- 
tion scale k, low-momentum fields are of the form (T^j)"*^ -1 ) (x) = 



Sj=o nfe'=j+i^A fc/ ^( x )- Successive i-derivations of scale k — 1, 
k — 2, . . . "push" (Tipi)^^ -1 ^ (x) downward like a down-going eleva- 
tor, in the sense that dt k _ 1 (Tipi)^^^ 1 ^ (x) is by construction of scale 

^x 

< k— 2, dt k _ 2 (T^j) - ^ -1 )^) of scale < k— 3 and so on. But of course, 

^x 

t-derivations may act on other fields instead. The last t-derivation act- 
ing on (T^j)^ k ~^> (x) drops (Ttpi)^^"^ (x) at a certain scale j < k. 
Then (Tipi)^^^ 1 ^ (x) leaves the elevator and is torn apart into its scale 
components ((T _> (- ? ' _1 ^i)' 1 ) h< . through free falling. Thus ip^ has a 
production scale k and a dropping scale j, while h itself may be called 
its free falling scale or simply its scale. 

(ii) If Pi > —D/2, then, at the dropping scale, (T^i) - ^ -1 ^) is separated 
from its average {Til)i)^^~ l \/\ x ) which must be dominated apart, 
while {Til)i)^(i~ l \x) — (T , 0j) -> (- J ~ 1 -'(Ai) splits into its scale compo- 
nents ((T^^ -1 ) 5 j i[)i) h ) h< - through free falling as in case (i). 
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The extension of the above procedure to lower scales is straightforward 
and leads to the following result. 

Definition 3.9 (multi-scale cluster expansion) 1. Fix a multi-scale 
horizontal cluster expansion = (F- 7 , . . . ,¥ p ) and consider a poly- 
mer down to scale j, P J ^, with horizontal skeleton F^. To such a 
polymer is associated a sum of products (G -monomial) v, (dressed in- 
teraction Ci n t{ ■ ]t)), where all t/^-variables such that A £ and 
ta < N ex t^ m ax + 0(n(A)) have been set to 0, and G is one of the mono- 



mials obtained by expanding 



f v dnt{ ■ \t)(x)dx 



n fc >,Vert fe .Hor* 

where Hor fc is as in Proposition \3.6\ and Vert fc is the following opera- 
tor, with N' ext max = N ext ,max + 0(n(A)), 

\r 4-k TT I aril f & ( 1_ *A) ext ' max p,Kxt, ma x 
Vert = || ^ d ^\t A =o+ dt ^^T- / —dt. 

a \ ^. — n •'0 



AeF fc \ r A =0 



(N' , - 1)\ tA 

\ ext, max ) 

(3.17) 



36 



2. (definition of n{A)) Fix and let A G W^. Then n(A) is the 
number of intervals of scale j(A) connected to A by the forest ¥^ A \ 

3. (definition of N(A)) Fix F- 7 ^ and some G -monomial, and let A G 
D J_> . Then N(A) is the number of fields tp^ A \x), x G A of scale 
J (A) /ymg in the interval A. 

Classically, one sets apart polymers made up of one isolated interval 
A 3 where no vertex has been produced; this means that t^j+i = for all 
intervals A^' +1 C A-?'; t AJ = 0; and s a ,- jA i = if A' G IP \ {A J }. Write as 
usual £j n t := X)g=i K q \ Kiq C q . Their contribution to the partition functions 
reads simply 

y d(i(^)e-f± c ^( ■ ;*=o)(*)«k 

= 1 — KqKq\ Kq J dfi(ip) ■ dv J C^ P { . ;t = 0)(x)dx • 
=: 1 - A 

by using a Taylor expansion to order 1. In the above expression, all fields 
have same scale j since all t-coefficients connecting A J have been set to zero. 
In all cases where the exponentiated interaction is bounded by 1 - which 
is the case in our (</>, d(p, <r)-model - the Cauchy-Schwarz inequality implies 
that \A | 2 < KX 2k , where k = min{K g ;g = l,...,p} since the integrated 
vertex is homogeneous of degree 0. Now, one eliminates the factor 1 by 
releasing the constraint that the disjoint union of all polymers must span 
the whole volume V, and ends up with a polymer whose evaluation is of 
order 0(X K ). This makes it possible to treat such polymers on an equal 
footing with the other polymers where some vertex has been produced (see 
introduction to section 6.1 for the general principles of the bounds for cluster 
expansions). 

Remark. Note that N(A) is at most of order 0(n(A)) for a single- 
scale cluster expansion. This is not true for a multi-scale cluster expansion. 
Namely, fix A h G B> h . Assume a low-momentum field ipi{x) (ft < —D/2) 
or 8^ip^(x) (ft > —D/2) has been produced at scale k and dropped inside 
A J at scale j, with k > j > h. Although the number of fields produced at 
scale k in an interval A h C A J increases exponentially with k — j, there are 
< N ex t,max + 0(n(A- J )) low-momentum fields dropped inside A J ' for a given 
Cx-monomial, since <?t . occurs at a power < N ex t )ma x + 0{n{A 3 )). On the 
other hand, the number of low- momentum fields ip h (x), x G A-?' with A-? C 
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A h originating from vertices of scale j may be of order #{A- ? 6 B J ;A J C 
A h } = M D ^- h \ This is a well-known phenomenon, called accumulation 
of low-momentum fields. This "negative" spring- factor must be combined 
with the rescaling spring factor M 2 ^- h \ see Corollary 12.71 resulting in 
j^(D+2/3)(j-h) ^ a p 0S m ve spring-factor if j3 < —D/2. This accounts for the 
(already mentioned) fact that secondary fields need not be produced for 
fields with scaling dimension < —D/2, see subsection 6.1 for details. 

At a given scale j, composing the horizontal cluster and momentum- 
decoupling expansions at all scales > j yields the following result, easy to 
show by induction: 

Lemma 3.10 (result of the expansion above scale j) 1. 

Z£(A;r^-i))= /d M (^- 1 ))^(A;(^) ft < J -_ 1 ;t^- 1 )) ) (3.19) 



with 



^r p (-H e 

\F3 EJFJ 



E 



(w p ) ( X £P ' X £P ) 



Yl Vert* Ho 

k>j 



-fy C(.;t)(x)dx 



(3.20) 



wh 



ere dfJ^Jt/j^) is a short-hand for Y\ P k=j d^ s k^ w k- } (ip k ). 

2. The right-hand side \3. 20\) depends on the low-momentum fields {^)h<j 
only through the dressed fields (Tij)i)~*^~ 1 \ since the t-variables of 
scale < j — 1 have not been touched. Hence it makes sense to consider 
the quantity Z'i^ p {\) := z(^ p (X; (^ fc ) fc <j_i = 0; i^- 1 ) = 1). 

3. The partition function Zy P (X) writes 



%? p w = E^T E Il^CPn) 

N ~ 1 non-overlapping Pi,...,P Jv e^~ > n ~ l 
oo N 

= e ^ e • n 



AT=1 



L,...,FjveP^ 71=1 



non— overlapping 

(3.21) 
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where Fhv> called polymer functional associated to horizontal (H) and 
vertical (V) cluster expansion, is the contribution of each polymer F ^ 
to the right-hand side of A3.20\) . 

Later on, the polymer functional Fhv wu 1 be replaced with the renor- 
malized (R) polymer functional Fhvr, or simply F, hence eq. (|3.2ip shall 
not be used in this form . The general discussion in the next subsection, 
which does not depend on the precise form of F, shows how to get rid of 
the non-overlapping conditions. 

3.4 Mayer expansion 

Recall the polynomial evaluation function Fjjv introduced in the expansion 
of Zy^ p , see eq. (|3.2ip . It is unsatisfactory in this form because of the 
non-overlapping conditions which make it impossible to compute directly a 
finite quantity out of it, see Introduction. Were this the only source of trou- 
ble, it would suffice to make a global Mayer expansion for Z v ~* p (\) after the 
multi-scale cluster expansion has been completed. However, it is also unsat- 
isfactory when renormalization is required; local parts of diverging graphs 
(for reasons accounted for in section 4) must be discarded and resummed 
into an exponential, thus leading to a counterterm in the interaction. This 
forces upon us a sequence of three moves at each scale j, starting from high- 
est scale p: (1) a separation of the local part of diverging graphs; (2) the 
Mayer expansion proper, at scale j; (3) the construction of the interaction 
counterterm (also called renormalization phase). 

Let us formalize this into the following: 

Induction hypothesis at scale j. After completing all expansions 
of scale > j + 1 and the horizontal/vertical cluster expansion at scale j, 
Zy P (X; t) has been rewritten as 

VCAji-^H f[ e |y|Mfc/fc " P(A) • [ dfi(^-V)Z v ^(\;t^-V;^ 
k=j+i J 

(3.22) 

where f k ~* p {\) may be reinterpreted as a scale k free energy density per 
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degree of freedom, with 

J2w 2 / dwJ " / d ^( w) (^) II F Hv(v n ; ip), 

N=l ' non-jr' -overlapping F 1 ,...,F N eV-'^ n=l 

(3.23) 

see Lemma \3.10\ for notations, where "non—j — overlapping" means thatWid 
W , . . . , P^v H D 3 are non-overlapping single-scale polymers, and Fjj V (¥ n ; tp) 
depends only on the values of ip on the support of P n . 

For j = p, this formula is equivalent to the single-scale expansion of 
Z^ p {\) of eq. EEZQ . 

According to the general expansion scheme (see introduction to section 
3), we must now perform three tasks. 

1 . Separation of local part of diverging graphs 

Consider the polymer evaluation function F HV {¥ n ; ip). If P n has N ext = 
^AieP„ninw r Ai < N ext ^ max external legs, situated in (possibly coin- 
ciding) intervals A], it may be superficially divergent, in which case 
one separates its local part according to the rule explained in section 
4: letting (tp^^ (%i))iel q i with xi G Af, be its external structure, 
F 3 HV {F n ;^) writes 

II / , ^4v,a m putated(Pn; *\ \{ 

(3.24) 

where 

4v,local(Pn^) = n //^kamputated^n^;^))- ( TF7 E ^ 

ie/, ^ V q] iei q 

(3.25) 

and F J HV (F n ; ip) minus its local part is now thought as if it had N ext; max — 
N ex t supplementary external legs shared in an arbitrary way among the 
intervals Aj, which produces a polymer denoted by S Next ' max ~ Next ¥ n 
belonging to V£ . 
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Equivalently (see subsection 4.1), considering only local parts, accord- 
ing to our new convention, 

4v(P»;0 = II /,^ F W,ampu t a t ed(IPn^;(x,)) " I[ 

l£lq » l£lq 

(3.26) 

where now 

^amputated (P-V'; 0*)) = TH E / II ^'^vO^ 0*))- 

(3.27) 

2. Mayer expansion at scale j 

We shall now apply the restricted cluster expansion, see Proposi- 
tion EH to the functional Z^ p (X;t; i/r+U- 1 )). The "objects" are 
multi-scale polymers P in O = {Pi, . . . ,P_/v} (see induction hypoth- 
esis above); a link t 6 L(0) is a pair of polymers {P n ,P ra /}, n / n'. 
Objects of type 2 are polymers with > N ex ^ ma x external legs whose 
non-overlap conditions may not be removed properly at this stage, 
because they are in any case dependent of each other in an as yet 
unforeseeable way through the low-momentum fields. All other ob- 
jects are of type 1, they belong to T >3 < ^r extimax ■= \& N e ^ max ~ Vj^; they 
are either vacuum polymers, i.e. polymers with no external legs, or 
superficially divergent polymers whose contribution to the interaction 
counterterm one would like to compute. 

Link weakenings S = (S {¥r} ) {¥r}eL{0) G [0, 1] L (°) act on Z^ P (X; t^O'" 1 ); ^rH?'- 1 )) 
by replacing the non-overlapping condition 

NonOverlap(Pl,...,Pjv) := ] J lp n ,P„, non-j -overlapping 

(Pn,P n /) 

= n n (i+(wa'-i)) 

(Pn,P„/) AeP„nDJ,A'eP„,nBJ 

(3.28) 

with a weakened non-overlapping condition 
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n n ^ 

(P n ,P„/) AeA e:Et (P„),A'eA e:rt (P n ,) 

1 + n ia ^ a 

\A6P„n]ttJ,A'eP n ,nDJ,(A,A')0Ae iC t(P„)xA e!ct (P n ,) 

where A ea; f (P) C P Pi W is the subset of intervals A J with external 
legs, i.e. such that r^j / 0. 

Note that each factor in (|3.29p ranges in [0, 1]. We ask the reader to 
accept this definition as it is and wait till the remark after Proposition 
13.111 for explanations. 

Let us now give some necessary precisions. The Mayer expansion 
is really applied to the non-overlap function NonOverlap and not to 
Zy^ p (\;t~ > ^~ 1 ';ip'^^~ 1 ^). Hence one must still extend the function 

/ dwi~> J V s( ^)(V^) Uti F J Hv( P n; 1>) to the case when the P n , n = 
1, . . . , N have some overlap. The natural way to do this is to assume 
that the random variables (V^ jp ) n =i....,N remain independent even 
when they overlap. This may be understood in the following way. 
Choose a different color for each polymer P n = Pi, . . . ,Pjv, and paint 
with that color all intervals A j £ P„fl W. If A j G A ext (P n ), then 
its external links to the interval A- 7-1 below it are left in black. The 
rules (|3,29p imply that intervals with different colors may superpose; 
on the other hand, external inclusion links may not, so that: (i) low- 
momentum fields t]j~*^~^{x), x € A J with A J G A ex t(P n ), do not 
superpose and may be left in black (till the next expansion stage at 
scale j — 1 at least); (ii) the color of high-momentum fields of scale j 
created at a later stage may be determined without ambiguity. 

Hence one must see ip 3 as living on a two-dimensional set, W x {colors}, 
so that copies of ip 3 with different colors are independent of each other. 
This defines a new, extended Gaussian measure dflgj^j^ip 3 ) associated 
to an extended field ipi : M. D x {colors} — > R, and Mayer- extended 
polymers. By abuse of notation, we shall skip the tilde in the sequel, 
and always implicitly extend the fields and the measures of each scale. 
Mayer-extended polymers shall be considered as (colored) polymers in 
section 6. 
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This gives the following expansion for Zy^ p (X; i - *^ 1 ^;tp^^ 

Proposition 3.11 (Mayer expansion) Let F{P^) be the set of all 

forests of polymers whose each component T is (i) either a tree of 
polymers of type 1 (called: unrooted treej; (ii) or a rooted tree of 
polymers such that only the root is of type 2. Then 



(3.30) 



with 

Mayer {Z{^ P {X; r^'" 1 ); ^O'" 1 )); F) 




(3.31) 

where Sg(W) is either or the minimum of the W-variables running 
along the unique path in ¥ from 0£ to o' v and F is the forest obtained 
from ¥ by merging all roots in max into a single vertex. 

As a result, (i) polymers Fg, F' £ linked by a Mayer link are j-over lapping 
(otherwise the derivative ds i would produce a zero factor) ; (ii) pairs of 
vacuum polymers F, F' belonging to different Mayer trees come with 
the factor 1: they have lost their non-overlap conditions and may su- 
perpose each other freely (in other words, they have become transpar- 
ent to each other). Hence Zy~* p (\; t - *^' -1 ); tp) factorizes as a product 



oo N 

Em E U fJ hvmK^) 

N=l non-j -overlapping ¥[,..., ¥' n eV^ n=l 



(3.32) 

where Fjj^ p M is the contribution of all (unrooted) trees of vacuum 
polymers. Denote by f J ~ > ' p (X) the quantity obtained by fixing one 
interval A J of scale j belonging to one of the polymers of the tree. 
Summing over all A J , one obtains an overall factor ^ V \ M1 f 3 ^ p W . 
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3. Renormalization phase 

For simplicity we shall assume that only 2-point functions (ifiiifii'}, 
1 < i, i' < d need to be renormalized (which is the case for the 
(0, dcf), a)-model). Denote by — \ f Aj dxbj i ,tp^^~ 1 \x)ip^^~ 1 \x) the 

contribution to Fj^ p M of all unrooted trees of polymers containing ex- 
actly one polymer with 2 external legs in A J , plus a cloud of vacuum 
polymers, attached to it, directly or indirectly, by Mayer links. The 
intervals without external legs of the different unrooted trees of poly- 
mers have become transparent to each other, and to the rooted trees 
too, hence b^, may be computed by considering one unrooted tree 
of polymers with 2 external legs, irrespectively of the position of the 
other trees of polymers. By translation invariance, fr^, is a constant, 
which is fixed by induction on j by demanding that 

(3.33) 

or better (so as to eliminate higher scale polymers which depend only 







(3.34) 

Hence eq. f|3.34j) yields \P iV as a sum over polymers (containing at 

least one interval of scale j) of an expression depending itself on IP- i / 
- an implicit equation. 

Local parts of 2-point functions are generated by the exponential over 

the whole volume, e 2 v *>* * ! , which may be seen 

as a counterterm in the interaction. This counterterm disappears by 
renormalizing El the Gaussian covariance kernel Cv, = Kq 1 - previ- 
ously renormalized down to scale j + 1 - to K , where 



(3.35) 



v 



13 Note that K - contrary to the original bare kernel - is only almost-diagonal, namely, 
J ip?(x)^i (x)dx — by momentum conservation if |supp(x J ) PI supp(x J )| = 0, i.e. if 

b'-j"l>2- 
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Note that the renormalization of C a in the case of the (0, deft, a) shall 
only be performed at scale p, while the other counterterms shall be 
left out as a counterterm in the interaction (see section 5). 

After applying a horizontal/ vertical expansion at scale to Zv~ 1 ) _>/3 ', 

we are finally back to the induction hypothesis, at scale (j — 1) this 
time. 

Let us now show the following single scale Mayer bounds. 

Proposition 3.12 (Mayer bounds) 1. (vacuum polymers ) Let Vq~*( A?) 
be the set of vacuum polymers down to scale j containing some fixed 
interval A J 6 W . Assume 

el p l- 1 |F^ y (P)|<K'A, (3.36) 

whereH F J HV (P) = f dp{il} j ^)F j HV (P; ^^). Then 

\f^(X)\<K'X(l + 0(X)) (3.37) 
with the same constant K' . 

2. (counterterm) Let V\i, (A 3 ) be the set of polymers down to scale j 

with exactly two external legs, ip~* ^ and ^\ displaced into 

the same point in a fixed interval A J . Rewrite F 3 HV tated {¥) := 

JM^fIwwmP-,^), p e Vi,; { An as 

^ i(1+A+ ^ ) i ; gS!^(A,P;^- 1 >), (3.38) 

with: 

• SpeP?'^(AJ),#{vertices of P}=2 F HV, amputated^ P ) = -^A 2 (1+0(A)); 

e r HV,amputated\ A i r ) 

PePf"^(Aj),#{vertices of P}>3 

= -K'X 3 {1 + O(A) + 0(M-^ 1+pi+p i'^,)). 

(3.39) 



14 Since P is a vacuum polymer, no high-momentum fields of scale j may be produced 
at a later stage, hence one may integrate out the field components ip 3 ^ . 
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Then 

6? ., = K\ 2 M j ( 1+ &+^(l + 0(A)). (3.40) 

We shall see in section 6.1 how to obtain estimates which are valid for a 
succession of Mayer expansions at each scale. 
Proof. 

1. Fix some interval A J £ W and compute / J ^ P (A) using Proposition 
EH] as 

EE E Mayer(Z^(A);T), 

N > 1 F 1 eV J Q ^ (A^),F 2 ,...,Pn&V ] ^ trees T over Plr-^JV 

(3.41) 

where 

1 * 

|Mayer(Z^(A);T)| < — ]J l^ry(^)- (3-42) 

' n=l 

The factor is matched by Cayley's theorem, which states that 
the number of trees over Pi , . . . , Pjv with fixed coordination number 
(?7.(Pi))j = i j ... i 7V equals pj (n(p-)-i)! • R- ecau that P^ and P^ are neces- 
sarily overlapping if £ 6 £(T). Start from the leaves and go down 
the branches of the tree inductively. Let Pi, . . . , P n (p')-i be the leaves 
attached onto one and the same vertex P' of T. Choose n(P') — 1 
(possibly non-distinct) vertices of P' n W (there are |P' n B J '| n ( p 
possibilities), fix their spatial location, A|, . . . , A^ij, and assume 

that A| G Pj PI W . For each choice of polymer P', this gives a supple- 
mentary factor < (if'A|P' D A j \) n ^ F to be multiplied by ^rerpm 
coming from Cayley's theorem. Summing over n(P') = 2, 3, . . ., yields 
e K A ' IP nAJ ' — I, which is < 2' p I (K'X) for A small enough. By induction, 
one gets |/- J ~ >P (A)| < J2h>i( K ' x ) h = K 'K l + 0{\)), where h is the 
height of the tree. 

2. The definition of 6^, and arguments analogous to those used in (1) 
yield (letting V^, := \~ 2 M~i ) tf.^, ) 

y>., = -if (1 + 0(A) + 0(X%,)), (3.43) 
hence the result by the implicit function theorem. 
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□ 

Remark. Consider the case when the Mayer expansion has surrounded 
P G V 3 ~* , a polymer with > N ex t >max external legs, by a cloud of vacuum 
polymers. Then the contribution of the cloud to Fj I y(¥;ip) may be com- 
puted as in Proposition l3.12l (1). see eq. f|3.37j) . leading at most to an overall 
mutliplicative coefficient < 1 + K'X(1 + O(X)) + . . . + (X ' A(1 +f (A))) " + . . . = 
1 + O(A) per interval A 6 P instead of e as in eq. (|3.36|) . This is important 
if A G A ex ^(P) is crossed by downward dashed lines as in §3.3 but contains 
no field, for an arbitrary product of factors e down to the scale of the low- 
momentum fields, possibly produced by successive Mayer expansions, may 
lead to a disaster. On the other hand (see end of §6.1.2) a factor 1 + O(A) 
per interval of P is not a problem. 

4 Power-counting and renormalization 

This section is divided into two parts. The first one gives a quick overview of 
how divergences in quantum field theory are discarded by extracting the local 
part of diverging diagrams; these ideas come from classical power- counting 
arguments which are recalled here. The whole idea of renormalization is then 
to transfer the sum of these local parts to the interaction as a counterterm 
(see §3.4). The second one is an informal discussion of the domination 
problem of low-momentum field averages, in particular in the case of the 
(4>, d(f>, cr)-model). The full treatment of this problem is postponed to §6.2. 

4.1 Power-counting and diverging graphs 

Definition 4.1 (Feynman diagrams) Let ip = (tpi(x), . . . , ipd(x)) be a 
multiscale Gaussian field with covariance kernel C^, and Ci n t = Yjq=i Kq^^^iq 
be an interaction. Then: 

1. a Feynman diagram for this theory is a connected graph V (whose lines, 
resp. vertices are generically denoted by I, resp. z ) with 

(i) external vertices y = (yi, . . . , y n ) of type 1, . . . ,p; 

(ii) internal vertices x of type 1, . . . ,p; 

(Hi) internal lines i connecting zi to z'^ with double index (i£,ig). 
Since the lines do not have a prefered orientation, and in order 
to avoid confusion, one writes I ~ (ig, zf, i' fJ z' £ ) or indifferently 
(i' e , z'fjii, z£). For every vertex z of type q, one may order the 
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lines leaving or ending in z as £ Z; i ~ (i±, z; i[, z[), . . . t z . n { z ) ~ 

so that n ( z ) = 1 and (*!'■■■)*«) = J q if z is 
an internal vertex, and n{z) < q, (i\, . . . , i q ) C I q if z is external. 

(iv) external lines £ ~ yf, i'g, y'^), where y' e G M D is an external 
point not belonging to F; \I q \ — n(y) of them per external vertex 
y of type q. The total number of external lines is N ext (T) : = 

J2q J2y of type q(\h\ ~ • 

The evaluation A(F) of an amputated Feynman diagram is given by 

Ap)(y)= [U dx ( II C f (i e ,z e ;i' e ,z' e )\, (4.1) 
J * \£eL int (r) J 

where x ranges over all internal vertices, and Li nt (T) is the set of all 
internal lines. 

The evaluation A{F) of a full Feynman diagram is given by 

A{T){y')= jj{dy[ \{ C^i e ,y e ;i' e ,y' e )\ A(F)(y), (4.2) 

where L ext (T) is the set of all external lines. 

2. A multi-scale Feynman diagram (r; (j(£))) is obtained from T by choos- 
ing a scale j(£) for each (internal or external) line and splitting each 
vertex into these different scales, with the following constraint, 

height(r) := mm{j(£);£ G L int (T)} - m&x{j(£);£ G L ext (T)} > 0. 

(4-3) 

Note that height (r) measures the height of internal lines of T with 
respect to the external lines it is attached to. 

The evaluation A(T) of a multi-scale Feynman diagram is given by 



A(T;(j(£)) eeLint{r) )(y) 



I 



U dx \ II C 3 /\i e ,z e ;i' e ,z , e )), 
e&L int (r) 



(4.4) 

while 
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A(r ] (j(£)heL { r ) W) = jU d y[ II c$ w (<M«;*£.«<)j^(r;0'W))(w) 

(4.5) 



^GLexi(r) 



Cluster expansions produce single-scale (horizontal) Feynman trees, muZ- 
tiplied by some G '-polynomial and by the exponential of the interaction, up to 
the modification of the measure by the weakening coefficients s. In the final 
bounds, see section 6, one bounds the exponential by a constant and applies 
the Cauchy-Schwarz inequality to the rest, yielding a multi-scale Feynman 
diagram. The scale j(z) of a vertex z is then the scale at which the vertex 
has been produced. 

Let us from now on assume that the interaction is just renormalizable. 
This means that, for every q = 1, . . . ,p, J2iei q ft = ~~ or in other words, 
J ipi q (x)dx is homogeneous of degree 0. In that case, the following very 
simple power-counting rules hold. 

Proposition 4.2 (power-counting) 1. The power- counting of an am- 
putated Feynman diagram V is the product A w ( r ) := A D \\ z A~ D Yie^L t (r) ^ 
see Definitional} for notations, where A is some large, indefinite con- 
stant representing an ultra-violet cut-off. 

Then the degree of divergence co(T) ofT is equal to D + J2eeL ext (r) Pit- 

2. The power- counting of a full multi-scale Feynman diagram T is the 
product 

2 teL(T) 
The multi-scale degree of divergence (jJ m . s .{T) is equal to 

oo m .sXT) = D • height(r) + E Mi*) ~ M)), (4-7) 

where L z is the set of internal or external lines leaving or ending in z. 

Let us give brief explanations. The power-counting in 1. may be obtained 
from Definition 12.41 by assuming all internal lines of T to be of scale 
and summing over all vertices except one, due to overall (approximate or 
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exact) translation invariance, see [39]. The multi-scale power counting is 
a refined version of the previous one, which takes into account the scale 
of the external legs; the definition of the height of T is somewhat ad-hoc 
and measures the horizontal freedom of movement of T with respect to its 
external legs if these are supposed to be fixed. More precise computations 
in the spirit of cluster expansion appear in the final bounds in section 6. 

The principle of the power-counting rule explained in section §6.1.2 is to 
rescale all fields produced at scale j of a multiscale Feynman diagram as if 
they were of scale j. The production scale of a given vertex z being j(z), 
this leads to a rescaled degree of divergence, 

<T led (r):=^.,.(r)-E E M{z)-m, (4.8) 

where the sum on I ranges over all internal lines leaving or ending in x. 
Hence 

< s f led (r) = ^-height(r) + ^ ^ pMy) _ m)} (49) 

where the sum on I ranges over all external lines leaving or ending in y. 

If height(r) and all j(y) — j(£) in (|4.9h are equal to the same positive 
constant, then o;^ s ^ aled (r) is proportional to the naive degree of divergence 
defined in Proposition 14.21 (1). 

Definition 4.3 (renormalization) // the degree of divergence w(r) of a 
Feynman diagram, resp. multiscale Feynman diagram, is > 0, then it needs 
to be renormalized. 

The local part of an amputated Feynman diagram or multiscale Feynman 
diagram is obtained by integrating over all external vertices except one (in 
order to take into account the global invariance by translation), 

A(T; (j(£)))(y = (2/1, • • • , Vn)) - j dy 2 ... dy n A(T; (i(£)))(y) =: Local(A(r ; (j (£)))). 

(4.10) 

By invariance by translation again, it does not depend on y\. 

Integrating over external points, one obtains the local part of a full Feyn- 
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man diagram, 



A(T; m))){y') - Local(i(r ; ■= 

Local(^(r ; (j(£)))) [Hdyl ]J cf\i e , y e : /, . ) 

' (4.11) 

which may be rewritten as 



n=[U d y( n s 

y V6i et( (r) 



Local(i(r;(jW))(y / )= l\\dy\ ff <%'>(i t , yr, i' e , y' e ) j A(T; (j0))(v)- 

(4.12) 

In oi/ier terms, all external legs ofT have been displaced at the same arbitrary 
external vertex, here y\. 

The renormalized amplitude 1ZA or 1ZA of a Feynman diagram or multi- 
scale Feynman diagram is the difference between its amplitude and its local 
part. 

Remark. Using a Fourier transform, the local part is equivalent to the 
classical evaluation at zero external momenta. 
A simple Taylor expansion of order 1 yields 

KA(T;(j(£)))(y') = jU d y E ( II CflfawM) 



y eeL ext (r) U'eL e xt(r),e'^£ 



L 



1 dt(y e - yi )E \ yi + t(y e - yi))^{y' t ) • .1(1: (jU ))(y). 



(4.13) 



In principle, this formula shows that the renormalized diagram amplitude 
1ZA(T; (j(£))) comes with a supplementary spring factor M _height ( r ), leading 
to an equivalent degree of divergence oj*(T) := oj(T) — 1. Namely, yi — y\ 
should be at most of order M~ mm ^^' ieLmt ^\ while by definition (see 



section 2) E [d^f\ Vl + t(y e - yiW^Hv'e) 
If uj(T) < 1 — D, then uj*(T) < —D and the diagram has become convergent 

im 



is of order AfM.M 



15 Otherwise one should renormalize to a higher order by removing the beginning of the 
Taylor expansion of the diagram in the external momenta. We shall not describe this 
straightforward extension of the procedure since we shall not use it in our context. 
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Let us give a formal proof of this general statement. This requires a 
little care because the covariance is not exponentially decreasing at large 
distances in our setting, so using systematically the Taylor expansion (|4.13p 
does not work. Consider a given multi-scale diagram (T; (jg)). We consider 
only divergent two-point subdiagrams for the sake of notations. Start from 
the divergent diagram of lowest scale, j m i n , say, with external legs y,y', 
and use the Taylor expansion (|4.13p . Choose a tree of propagators of scale 
> j min connecting y to y' through intermediary points X2, ■ ■ ■ ,x n , and set 
x\ = y,x n+ i = y', so that \y - y'\ < Yl m =i \x m - x m+1 \. The divergent 
diagram with external legs y, y' has been renormalized as explained above. 
Letting j m be the scale of the link (x m ,x m+ i), one has lost a rescaled fac- 
tor J2m=i djmm (xm,x m+1 ) = J2 m =i M ~ {jm ~ jmm) d jm (x m ,x m+1 ). Each term 
in this sum has obtained a spring factor M~^ m ~^ min \ which shows that 
the corresponding diag ram with, external legs is already de facto 

superficially renormalized, although subdiagrams of higher scale may still 
need renormalization. In the process, it is clear that every propagator of the 
diagram may appear only in one tree of propagators at most. In the bounds 
of section 6, this yields an overall factor per polymer P which is bounded 
W rbei(P)(l + d ( x ti x'p)), which is easily controlled by the polynomial 
decrease of the covariance at large distances. 

Let us summarize our brief discussion. 

Definition 4.4 (diverging graphs) A Feynman graph or multi-scale Feyn- 
man graph T is divergent if and only if 

uj(r):=D+ Y, Pk>°- ( 4 - 14 ) 

We shall call N ex t tma x the minimum value of N ex t such that every dia- 
gram with > N ex t^ m ax external legs is convergent 

Renormalizing the evaluation of a Feynman graph or multi-scale Feyn- 
man graph V yields a quantity 1ZA{T) for which the displaced external legs 
have obtained a supplementary spring factor, which is equivalent to replacing 
one of the scaling dimensions /3i t , £ € L ext (T) with (3* := j3i t — 1, or globally 
w(r) by uj*(T) :=uj(r) - 1. 

Let oj max < the maximal value of the set {oj*(T)}, where T ranges over 
the set of all Feynman graphs. 

Example ((^>, d(f>, cr)-model) 

In our case (3$ = a, Pq^ = a — 1, (3 a = —2a. Note however that 
the interaction (see section 5) is such that split vertices are either of type 
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d(ftj(j) k cr k or o~i d(fi k 4> k , with k' = k ± 1, hence (assuming k' = k which does 
not change anything) <j) may not be an external leg of a multi-scale Feynman 
diagram (in particular, the vertex (j)d(pa is not renormalized). Let Ng^, resp. 
N a be the number of external d(p- or c-legs of such a diagram T. The power- 
counting rule yields YlieL ex t(T) Pit = ( a ~ ^)^dcj> ~ 2azN a < — 1 as soon as 
Nd<j> > 2 or A^ > 4, so N ext: max = 4. However, by symmetry arguments, 
local parts of diagrams with Nq^ or N a odd vanish, so there remains only 
the cr-propagator (<7±(x)cr±(y)) to renormalize. 

Let us write down precisely these power-counting rules in the framework 
of multi-scale cluster expansion, where Feynman diagrams are replaced by 
polymers. Recall from §3.3 that a low momentum field (J 1- ^ -1 )^)' 1 ^) - 
or {T^^^ 5^i) h {x) - has three scales attached to it. It is the difference 
between the two highest ones - the production scale k and the dropping scale 
j, with k > j - that counts here. Namely, considering h < j' < j , there are 
two cases: 

(i) either t . j> (the number of i-derivatives produced inside A x ) is > 
N e xt,max- Then rescaling the fields to which the i-derivatives are ap- 
plied is enough to ensure the convergence of the polymer W ~* con- 
taining A x ; 

(ii) or t j> < N ex t^ m ax- Then (by definition) t j/ is set to 0, which kills 

l^x ^x 

{T^^~^'^)i) h {x) or (T^^~^S^il)i) h (x). Hence one must count only on 
the derived fields for the power-counting. As explained above, this 
may give a non-negative degree of divergence, in which case one must 
renormalize by removing the local part of the polymer. 

This means that the rescaling spring factor M^^ k ~ h ^ of these low-momentum 
fields must be split into M^ k ~^ - ensuring the horizontal fixing of the poly- 
mer, and counting in the right-hand side of eq. f^.ffp - and M^^~ h ' - which 
helps control the accumulation of low-momentum fields as explained briefly 
in 13.3. 

4.2 Domination problem and boundary term in the interac- 
tion 

Assume ipi is a multi-scale Gaussian field with scaling dimension G 
(— D/2,0). Then 2/3j > —D so, by Definition I4.4| the associated two- 
point functions must be renormalized. If tpi occurs in some interaction term 
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X Kq ipi q , i.e. i € I q , then renormalization produces a scale-dependent coun- 
terterm of the form 

SC^i-x) = J2^f q C J M^ D+2 ^ (x), (4.15) 

where \j is the running coupling constant, and Cj is a scale-dependent con- 
stant. In this paper we generally assume that Xj = A is not renormalized. 
(In the case of our (<fi, d<j), <r)-model, SC is of the form Sj>o ^(O^ 7 ) - ^) 2 ' 
with hP ~ X 2 M^ 1 ~ 4:a ^ and Cj ~ 1. See the counterterm 5C4 in section 5 be- 
low for details.) Separating the low-momentum field ^(x), x S A J into 

(x) as in Definition 12.51 produces a secondary field 

S^ipJ^ ^ whose contributions to the partition function are easily bounded 
thanks to the "spring factor" (see subsection 6.1). Unfortunately, the av- 
eraged fields ipj*^ (A J ) do not come with such a spring factor and (un- 
less /3i < —D/2, see subsection 6.1 again) must be dominated apart, using 
the positivity of the interaction. We assume that Cj > 0. Then Lemma 
E32 (1) shows that |A K s (T^)^ '^ (A-»')| n ( A ) e~ / ' S£-{^;x)dx ig bounded by 

/ _i \n(A) 

O I i^A)^- 2 , which agrees - up to unessential local facto- 

rials n(A) n ( A )/ 2 , see §6.1 - with the correct power-counting, but the "petit 
facteur" \ K - q has been entirely used, in contradiction with the general guide- 

line of cluster expansions. So - unless C- 2 is small - one must use a different 
strategy. 

Several strategies have been used, depending on the model. Here things 
are particularly simple because the mass counterterm of highest scale, b p ~ 
\ 2 M p<yl ~ Aa \ couples with all scales of the field a and is much better than 
the term of order A 2 M^ 1_4a - ) appearing in the right-hand side of 14.151 This 
is the reason why this counterterm has been set apart from the counterterm 
and put into the covariance of a right from the beginning. The supple- 
mentary spring factor M - ^ 1-4 ")^ - -?) per field plays the role of a "petit 
facteur", except in a small scale interval p — q,. . . , p, where q ~ ln(l/A) 
is p-independent. Add now e.g. to the interaction a term of the form 
M -(4na-i) PX KY ip , <p \\(Ta)^P'(x)\\ 2n with 2n > 4, which is homogeneous 
of degree 0, and totally negligible away from the highest scales because 
of the evanescent coupling coefficient M~( 4na_1 ) p < M~^ a ^ p < 1. If 
k < 2n, then each c-field in the interaction is coupled to \ K / 2n ^> A, which 
makes it possible to dominate the low-momentum fields for the highest scales 
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p — q + 1, . . . , p. In this sense this term is a (Fourier) boundary term. The 
term 5L\2 in section 5 below plays this role. The choice of 4n = 12 is 
arbitrary. 

5 Definition of the model 

As a general rule, we shall denote by px the Gaussian measure with covari- 
ance K^ 1 if K is a positive-definite kernel. 

Recall from the heuristics in section 1.2 that one wants to define as 
interaction J Ci n t((p,a)(x)dx = iX J (dA + (x)a + (x) — dA~ (x)a-(x))dx. In- 
tegrating over R and splitting the fields (pi,4>2, °~± into their different mo- 
mentum scales yields a momentum conservation rule; as a result, products 
<9</>-f (x)4>2 '{x)a\(x) (j\ < 3i) or d^(x)^{x)u h _(x) (ji < ji) contribute to 
the action integral only if the two highest scales differ by at most 1 (pro- 
vided M > 2, say). The notation YT^j 2 k means that the sum is restricted to 

this admissible subset; explicitly, Z)(^j'2,fc)e/(' ' ' ) = ^(j u j2,k)eim adm > witn 
jadm = j 2; k); ji ~ J2>k or j 2 ~ k > ji or k ~ ji > j 2 }, where j ~ k 
means j = k or k ± 1. This scale restriction is added explicitly by hand 
because the dressing procedure (which is not translation-invariant) spoils 
momentum conservation. 

Definition 5.1 (propagators) (i) Let dpT rp {(t>) = dp,\^\i+2 t x{(f)^' p ) be the 
Gaussian measure associated to the field (p^ p defined in section 2.2. 

(ii) Let dpT* p {a±) = dii\^i-Aai A+hP {a± p ) be the Gaussian measure associ- 
ated to the fields a ± p defined in section 2.2. 

The two-by-two matrix coefficient b p is called the renormalized mass co- 
efficient of the ex-field at scale p. It is equal to the local part at scale p 
of the two-point function of the a- field (see precise definition below). Note 
that dfi^i-4 aM+bP (a^ p ) = -^re^2 bP S(^ p ) 2 i x ) dx d^i-ic, where Z' is a nor- 
malization constant. This quadratic term appears with the opposite sign in 
the interaction Lagrangian £~^ p ((j), o-)(x) below, so the interacting measure 
has really been constructed initially by using a <r-field with bare covariance 
d//|£|i-4a, as explained in section 1. 

Note that all summations over scales in the interaction start from the 
scale j = 1. Namely, it is really meant to cure ultra-violet divergence and 
absolutely useless in the infra-red region. However, the cr-field is not coupled 
to the infra-red part of the Levy area, which yields a regularized area which 
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is not homogeneous any more. This can easily be corrected by extending 
the Fourier partition (x^)j>o to the infra-red region |£| < 1, with |suppx J | C 
[M- 7-1 , AP +1 ] for all j £ Z. One may then sum over all scales from j = p to 
— p, without bothering to renormalize at scales j < 0, and let p — > oo. We 
leave details to the reader and concentrate on the region |£| > 1. 

In the sequel, expressions such as ba 2 are to be understood as a scalar 
product (ba,a) = 6 +i+ (cr + ) 2 + 26 +i _(j + <7_ + b--(a-) 2 . 

Definition 5.2 (interaction) (i) Let 

hP 

q%{M(x) := Cf'faaXx) - -(a^(x)) 2 , (5.1) 



where 

C^"{4>,a)(x) 



\ D +\ E 8^W4 2 W^W " D - ^ {x)4>{\x)v k _(x)\\ 

V \l<i 1 ,i 2 ,k<P J \l<j 1 ,j 2 ,k<P j) 

C adm adm ") 

i\i J2 d4>{ x (x)4, 2 2 (x)a k + (x) - ]T 84,^ (x)^ 1 (x)a k _(x) [(5.2) 

j\ f adm adm \ 

-I ]T d4>{(x)^x)a k + (x) - ]T a^W^W^fi) ; (5.3) 

|.l<j<p,l<fe<p l<j<p,l<fc<p J 



(m) (dressed interaction) 
Let 



• ;t)( X )--(t^\\(Tar^x)\\ 2 + 5C^( . ;t)(x) + 5C^( . ;t)(x), 

(5.4) 



where: 

Ct p ( • ;t)(x)=£jf( • ;t)(x)-C^( . ;t)(x), (5.5) 



£^:=UB+ a(T^ 1 )«(:r)(T^ 2 ) 3 2(:r)(T^< 3 ( T + )' e ( ;c ) 
2<p'<p Kil.B.'Sf'- 1 
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and similarly for C^ p (the Fourier projections D± have been defined 
in section 1); 

!*)(*):= 5 E hP '~ l E (l-{€) 2 )({Tar^"-^{x)) 2 - 

2<P'<P 2<p"<p' 

(5.7) 



;t)(s) := 

M -(12a-l)p A 3 J ||( rff )^P( x )||6 + V (l-(^') 6 ) / IKTct)^'- 1 )^)!) 6 ^ 
I 2<p'<p • /A 



(5.8) 



/or i/ie Euclidean norm \ \a(x)\\ = y/\a+(x)\ 2 + |o~_(x)| 2 . 
5y construction C^(<p,a) = £^(<p,a;t = 1). 

(m) Let 

Z^ p {\):= J e-fv^i^^i^d^P^d^Pia), (5.9) 
and, more generally, 

Z^ p (\;(0 h ) h<p ,,(o- h ) h<p ,) = J e~fv^(^t=m^ d ^^ P ^ )dfx P'^ P{a) 

(5.10) 

which is a function of the low-momentum components of the fields, 
considered as external fields. 

(iv) (definition of the renormalized mass coefficient b p> ) Fixb p ' by requiring 
that 

f d 2 i 

= J d y dtf+tf-i) Mdo-^iP 1 - 1 ) (y) Z v~* P ( X > & W' ( a W)L^(p'-i)=^c p '-i) =0 - 

(5.11) 

Note that 5C^ P ( . ; t)(x) may be seen as a dressed interaction as in Def- 
inition 13.81 if one sets b p := (the counterterm of scale p has been treated 
separately). The interaction before dressing therefore vanishes, which is 
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coherent with the fact that it has not been put into the model from the be- 
ginning, but built inductively to compensate local parts of diverging graphs. 

Although the actual form of £^ p looks complicated, it is very close to 
the interaction term iA (dA + (x)a+(x) — dA~(x)a-(x)) of section 1, see eq. 

dm, dm, {ttzd. 

6 Bounds 

6.1 Gaussian bounds 

This paragraph is the backbone of the section, since it provides a means 
(i) to bound the sum all possible Wick pairings of all possible G-monomial 
associated to a given polymer; (ii) to bound the sum over all possible poly- 
mers containing some fixed interval A at its lowest scale. The general idea 
(as explained in the introduction) is that a polymer is connected either by 
horizontal cluster links which are polynomially decreasing at large distances, 
or by vertical inclusion links which create spring factors. The computations 
below for a polymer P (see §6.1.2 and 6.1.3) give in the end a bound which is 
of order Y\ v X K YIa M~ £ for some positive exponents n, e, where the product 
ranges over all vertices v and intervals A of P. This is the general prin- 
ciple of the bounds for cluster expansions. Both terms X K and M~ £ are 
called a "petit facteur par carre" (small factor per interval, in French). It 
does not include combinatorial factors (see §6.1.4), dominated averaged low- 
momentum fields (see §6.2), and possibly some other terms (see first step in 
the Proof of Theorem 6.2 in §6.3), but all these are proved to give for each 
choice of cluster expansion and G-polynomial a supplementary factor of the 
type IIaU + 0(X K )) for some exponent k' > 0, which is eaten up by the 
factor J]a m ~ S - 

Since all computations are Gaussian in this paragraph, we shall take the 
liberty to write ( (• • • ) ) instead of E[ (• • • ) ], without any risk of confusion. 

6.1.1 Wick's formula and applications 

We first recall the classical Wick formula. 

Proposition 6.1 (Wick's formula) Let (X%, . . . ,X2n) be a (centered) Gaus- 
sian vector. Pair the indices 1, . . . , 2N; the result may be represented as a 
graph ¥ with n connected components, linking the vertices 1, . . . ,2N two by 
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two. As in Definition 13.51 we use the pair notation I = {ie,i' e } for links. 
Then 

(X 1 ...X 2N )= Yl H(XuXi' t )- (6- 1 ) 

F pairing of {1,...,2N} i& 



Proof, see e.g. [E], §5.1.2. 



□ 



Corollary 6.2 (simple Wick bound) Let (X%, . . . , X 2 n) be a Gaussian 
vector. Then, for every K > 0, 



27V-1 



\(x 1 ...x 2N )\<k- n n 



1=1 



j>i 



•2) 



In particular, 



2N 



\(xi...x 2N )\<n 



i=i 



l + J2\(X t Xj 



.3) 



Proof. Expand the right-hand side and use Wick's formula to get eq. 
(|6.2p with K = 1. The bound with K 7^ 1 may be obtained from the 
previous one by a simple rescaling X% — > \[KXi, % = 1, . . . , 2iV. □ 

The above bound eq. (|6.2p depends on the ordering of the variables 
Xi, . . . ,X 2 jy, although {X\ . . . X 2 n) doesn't, of course. The idea conveyed 
by this bound is that it may be important to choose the right order. Simi- 
larly, eq. ()6.2p is clearly optimal when the factors K^2j >i \ (XiXj)\, 1 < i < 
2N — 1, are of order 1. 

However this bound is too simple to apply in most cases, and we shall 
need refined versions of it using the spatial structure of the Gaussian vari- 
ables. The following lemmas, for the reasons we have just explained, are to 
be used after a suitable rescaling. 

Corollary 6.3 (Wick bound with spatial structure) 1. (single-scale 
bound) Let (-X'(A,n))AeD',l<n<iV m oa(A) be a Gaussian vector indexed by 
M-adic intervals of scale j. Denote by I = {(A,n);A £ 3 ,1 < 
n < iV(A)} the total set of indices. Call connecting pairing a partial 
pairing ¥ of the indices (A, n) such that its spatial projection F with 
vertices {A € W; 3n < N max (A) \ (A,n) € F} and links {{A, A'} G 
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W x W;3n < N max (A),n' < N max (A') | (A,n) ~ F (A',n')} is con- 
nected. Fix some M-adic interval Ai G B- 7 . T/ien 

£ ,n^(A^))l 

F connecting pairing of 7,|F|=2JV,AieF feF 

/ /N max (A) ^ ^ 3N 

" ( 1 + f U i? ( E E K*(A,n)*(A',n')>l 

n=l (A',n')e/,(A',n')^(A,n) 

(6.4) 



A ecu 



2. (single-scale bound, improved version) 

More generally, assume N max (A) = oo and K : (W x {1, . . . , d}) x 
(W x {1, . . . , d}) — > R+ is some kernel, copied an infinite number of 
times, so that K((A,pd + i), (A',p'd + i')) = K((A,i), (A',i')), with 
p,p' = 1, 2, . . . (see Remark below). Let I = W ' x {1, 2, . . .}. Conect- 
ing pairings of I must be understood modulo the pair identifications 
((A,pd + i), (A',p'd + i')) ~ ((A, »), (A', »'))■ 

TTien, /or any 7 > 1, letting N(A) := #{(A,n); (A,n) G F} - to 6e 
interpreted as the number of fields lying in a given interval A -, 

£ J] (N(A e )N(A' e ))^ K((A e , ie ), (A' e ,i' e )) 

F connecting pairing of /,|F|=2JV,AieF feF 

< [l + sup Yl E ^((A,i),(A',i')) 

V AeBJ A'GD^ i,i'=l 

(6.5) 

3. (multi- scale bound) 

Let K : (W^ x {1, . . . , d}) x (D-^ x {1, . . . , d}) ->• M+ 6e some ZcerneZ 
indexed by M-adic intervals of scale > j, copied an infinite number of 
times as in 2. We denote once again by L the total set of indices. Let 
L k := {(A,n) G /; A G B k }, k>j and L~> k ' := Wj' =j J-»", fc' > j. Fix a 
certain number of (non necessarily distinct) M-adic intervals for each 
scale k = j,j + l,...,p, say, A k ,..., A k k (k > j), with c k > (k > j) 
and Cj = 1; write for short A 3 ^ = {{A k )k>j,i< c <c k } C . Let 
F 3 ^(A 3 ^) be the set of multi-scale cluster forests F 3 ^ (called: A 3 ^- 
connected multiscale cluster forests,) such that, for each j' > j, each 
vertex ofF 3 is connected by horizontal cluster links or inclusion links 
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to some (possibly many) of the selected intervals (A k )k'>j',l<c<c k is an< ^ 

are not connected within F J by horizon- 



the intervals A 



i ) 



A J 



tal cluster nor inclusion links. (In other words, each selected interval 
Aj , . . . , Ac , lies within a different horizontal cluster and inclusion 

connected component o/F 3 ''"^, and these cy connected components ex- 
haust the set of horizontal cluster and inclusion connected components 
of F J ~* which contain at least one M-adic interval of scale f). Call 
A- 7 "^-connecting pairing a partial pairing F of the indices (A, n) such 
that its spatial projection F has the same set of vertices and links as 
some A- 5 -connected forest, plus possibly some supplementary links, 
possibly reducing the number of connected components. Then : 



H (JV(Aj)AT(AJ)) 7 tf((A,,i,),(A 



FA^ —connecting pairing of I,\¥\=2N?& 

d 

< | 1 + max sup 

k>j A fc eD fe 



E £ K((A\i),(A>,i>)) 

A'eW^ k i,i'=l 
fc-1 d 1 X 3JV 

E E E ^((aVmav)) 



where A k> is the unique interval of scale k' such that A k ' D A k . 

Remark. When N max (A) = oo, which is due to the fact that the 
total number of fields in a given M-adic interval, iV(A), may be of order 
n(A), hence unbounded, the bound in 1. is infinite. In practice, a cluster 
expansion generates - thanks to the polynomial decrease in the distance of 
the covariance of multi-scale Gaussian fields - extremely small factors per 
interval when n(A) is large. The idea is then to bound \ {ipl (x)i[r!,(x')}\, x £ 
A,x' £ A' by i), (A',i')), where some of the polynomial 

decrease in the distance has been retained in the kernel K (see §6.1.2). 

Proof. 

1 . Consider first the left-hand side of (|6.4p . Consider a connecting pairing 
F such that |F| = 2n and containing the M-adic interval Ai, and a 
spanning tree T of F containing A^. Associate to F the following 
sequence of links and of factors 1: 

- consider all the pairings of the indices (Ai,n), 1 < n < N max {A{) 
among themselves and with indices (A',n'), A' ^ A\; say (in some 
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arbitrary order), (Ai,ni) pairs with (A^n^), . . ., (Ai, njVi-i) pairs 
with (A^ Vi _ 1 , tijy J. Insert after these Aq — 1 links a factor 1, signi- 
fying that all paired Gaussian variables lying in the interval Ai have 
been exhausted; 

- continue to explore new vertices of F by going along the branches 
of T. Always insert a factor 1 after all the pairings of the Gaussian 
variables lying in a given interval have been exhausted. 

Since T is connected, all M-adic intervals in F and all indices in F will 
eventually have been explored. The number of factors is £(¥) + |F| < 
N + |F| = 3N, to the completed by the required number of factors 1 
so that there are exactly 3iV factors. 

Consider now the right-hand side. Let 

Ka:= £ £ |<*(A,„)*(A>0>l (6-7) 

l<n<N max (A) (A',ra')^(A,n) 

and K$ = 1, and expand K 3N := (Kq + sup AgD j K&) 3N . One gets 

Km = E E (-p^a)^- 1 -i • 

V Ni<...<N p A 

(sup^A)^-^" 1 • (sup^A)^-^- 1-1 • 

A A 

(6.8) 

Replace the first sequence of Aq — 1 factors by K^ 1 ^ 1 < (sup A K^) 1 ^ 1 ^ 1 
and expand them, which encodes in particular all possible pairings of 
the Gaussian variables lying in Ai. Consider now an interval A2 7^ Ai 
linked by T to Ai, and replace the second sequence of factors by 
K^~ Nl_1 , and so on. Thus one has encoded all possible connect- 
ing pairings containing the interval Ai. 

2. Considering as in the previous case all the pairings of the indices 
(Ai, n), n = 1, 2, ... for a given pairing F, one gets the factor 

K Al ■= £ ((N(A 1 )N(A[))^K((A 1 ,i 1 ), (A[,i[))) . . . 

((7V(A 1 )iV(A^ 1 _ 1 ))-^((A 1 ,i i v 1 _ 1 ),(A^ 1 _ 1 ,4 1 _ 1 ))), 

(6.9) 

where the sum ranges over all possible (Aq— l)-uple couplings ((A 1; i), (A' 
with multiplicities (note that < Aq — 1 < iV(Ai), depending on 

the number of couplings of fields inside Ai). Then 
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, iV(Ai) N(A') i Nl 1 

«ax < E E ^((Ai,i),(A',o) E E WAiWAOr 

\A'eFJ l<i,i'<d n=l n'=l 

< EE ^((Al 0,(^,0)) • (6.io) 

Apart from this slight difference, the exploration procedure is the 
same. 

3. Choose a spanning tree of FnB p , complete it into a spanning tree of Fn 
jjCp- 1 )-^ anc i go on . As in 1., explore the horizontal cluster connected 
components of scale p starting from the selected intervals Ac, 1 < c < 
Cp, then the connected components of scale p — 1, and so on, down 
to scale j. The only difference is that two different horizontal cluster 
connected components of F of the same scale f may be connected from 
above by inclusion links and horizontal cluster links of higher scale; in 
this case, this procedure may not explore all vertices of F. Fortunately, 
the bound in eq. (|6,6p gives the possibility, starting from some interval 
A fc g B fc , to go on to explore all the Gaussian variables located below 
A fc , i.e. in some interval A fc D A k with k' < k. 



□ 



6.1.2 Gaussian bounds for cluster expansions 

We assume here that Ci n t is just renormalizable, so that (assuming just for 
simplicity of notations that its coefficients are scale-independent) Ci n t = 
K\ A Kl + ... + K p X K Pifj Ip , where k u . . . , K p > and J2ieh Pi = ■ ■ ■ = 
^2i£i Pi = — D. Each term tp^ , . . . , ipi p is called a vertex by reference to the 
Feynman diagram representation (see section 4). Let us recall briefly that 
the G-monomials are produced: 

- either by horizontal cluster expansions; if i£ S I q , £ being a link at 
scale j, then ^ e~ XKq J^ I i^ dx produces \ Kq ip Iq \{ ie }(xe). On the other 

hand, — t-, — - may derivate the low-momentum components of monomials 

produced at scales > j + 1, which lowers the degree of G; 

- or by t-derivations acting on e~ J c -mt{ x ) dx ; yielding (up to t-coefficients) 
some (scale components) of the \ Kq ipi q (xA), 1 < Q < p, integrated over 
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some M-adic interval A. Again, t-derivations may derivate the monomials 
produced at scales > j, which does not change the degree of G. 

The above products of fields must now be split according to their scale 
decomposition. Thus one obtains a certain number v of vertices split into 
different scales. 

We suggest the following notations in order to avoid the proliferation of 
indices. Make a list (tpi, . . . ,ipd), d = \I\\ + . . . + \I P \, of all the fields involved 
in the interaction, possibly with repetitions. Thus the cluster expansion 
at scale j generates at the same scale X Kq ipj (xg), where Ig C I q \ {ig} C 
{1, . . . , d} \ {ig} for some q < p; on the other hand, each t-derivation in an 
interval A G W generates (up to i-coefficients) some X Kq if) J lA (x&), C I q , or 

(with an extra index ta for the order of derivation) (X Kq ip 3 lA (xa,t))t=i,...,t a • 
But other field components of scale j, lying in some fixed interval A J G 
W, are produced, either at an earlier stage k > j, in the form of a low- 
momentum field, ipi{x) or 8 k ^{x) (secondary field) with x £ A k , A k G D fe , 
A fe C A J , or at a later stage h < j, in the form of a high-momentum field, 
Res h A ,i; j {x), x G A h , where A h €B h , A h D A*. 

The general principle of bounds for cluster expansions in quantum field 
theory (as explained at the beginning of §6.1) is to (1) use the polynomial 
decrease in the distance of the covariance of the field components; (2) find 
out a "petit facteur par carre" (small factor per cube, or rather per interval 
in one dimension). This means essentially the following: chose some possibly 
derivated interaction term A K9 V/,\{^}(^)) x e G A J or X K i^j q (x A j) coming 
from a vertex at scale j; the fields i G I q scale like M~^ iJ , and the 
integration over the interval A J G B- 7 produces a factor M~i (or M~ D i in 
general). As for the cluster expansion at scale j, it has produced a factor 
C^(i£,xf,i' e ,x' e ) which scales like , times the same quantity with a 

prime. Supposing one chooses a splitting of the vertex such that all fields 
are of scale j, then the product of these factors is X Kq M~ Dj M ^^tq^ 1 = 
X Kq <C 1, which is the "petit facteur". Unfortunately the splittings of the 
vertex produce much more complicated situations; however, the guideline 
is to compare the scalings of the high-momentum fields (rewritten as a sum 
of restricted fields) and of the low-momentum fields (possibly rewritten as 
secondary fields, modulo averaged fields) with the scaling they would produce 
if they were of scale j. The high-, resp. low-momentum rescaled fields, write 

^{x) =: M-^ fe -^> rescaIed (x) (* > j), resp. ^(x) =: M^- h ^' ieBtxiBd (x) (h < j), 

(6.11) 
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, iq is accompanied 



yielding positive or negative rescaling spring-factors (depending also on the 
sign of /?), see after Corollary 12.71 The factor X Kq is split into the different 
scales of the vertex, so that each field ipi(x), i = ii, 
by a small factor < A K °/I /<3 I for some kq > 0. 

Averaged fields must be treated apart and account for the so-called dom- 
ination problem; bounding them may require part of the small factor X Kq , so 
that, generally speaking, the "petit facteur" is of order A K , for some k > 
but small. 

Let us first consider the following single scale situation, throwing away 
all low- or high-momentum fields for the time being. 

Lemma 6.4 Let tp = (tpi(x), . . . ,ipd(%)) be a Gaussian field with d compo- 
nents such that 



\CUi, X] i',x')\ = \(i4(x)Mx'))\<K r 



(1 + W\x - x'\ 



(6.12) 



for every r > 0, with some constant K r depending only on r; these bounds 
hold in particular if ipi, . . . ,ipd C {(<9 n, 0i)n>Oi ••• > {d n i>d)n>o} are deriva- 
tives of some independent multiscale Gaussian fields ipi, . . . 

Consider a horizontal cluster forest F 3 6 T 3 of scale j, and associated 
cluster points X£,x' e , £ £ L(F 3 ), X£ £ Ag, x' e G A' £ . Choose: 

- for each link i £ L(F 3 ), a subset I# of {1, . . . ,d}\ {ie} and a subset 1'^ 
of{l,...,d}\{i' e }; 

- for each M-adic interval A G F 3 , T/\ < Next,max 

+ 0(n(A)) subsets 

(-^A,t)t=1 ... ta ■ ■ ■ ) d>}, an d additional integration points (xa,t)t=1 ... ta 

in A. 

Such a choice defines uniquely a monomial G 3,3 = G 3 ' 3 (F 3 ; (Ie), (i^); (Ia,t), (xa,i 
in the fields ipi, i = 1, ... ,d taken at the cluster points (xg), (x'«) and the 
t-derivation points (xa,t)> namely, 



G> 



,3 .- 



yi Cj v{y ':(: :: ) 



^GL(Fi) 



II no* 

AgFJ t=1 



(6-13) 

where v(F 3 ;G 3 ' 3 ) := 2L(F 3 ) + X^AeF^ t a *s ^ e total number of vertices 
obtained from the horizontal cluster and the t-derivations is (two per hor- 
izontal cluster link, one per t-derivation acting on the exponential of the 
interaction). 

Denote by N 3 (G 3 ' 3 ; A),i = 1, . . . , d the number of fields ipj (x), x £ A oc- 
curring in G 3 ' 3 if A e F 3 , summing up to N 3 (G 3 > 3 ; A) := J2Li N H Gjd '> A )> 



65 



and by N- (G 3 ' J ) := X^Agfj ^1 A) the total number of fields ipj occur- 
ring in G hl . Similarly, we denote by N- (¥ J ) the number of half-propagators 
of type i in F 3 . Note that N j (G j ' j ;A) = 0(n(A)), see W13\) . 

1. Let 

^au Ssi an(^;^):= !d^) \\ C^Xf^x',) ■ G 3 ' 3 . (6.14) 

£eL(Wi) 

Then 

d 

\I j . (¥ j - G j ' j )\ < K^^X Kov( - ¥j ' G3 '^ TT M- j ^ N ^ G3 ' 3)+N ^ w)) . 

I GcLUSSiclIl v ' /I — JL JL 

i=l 

(6.15) 

2. Fix the total number of vertices, v = v(¥ 3 ;G 3 ' 3 ), and fix one M-adic 
interval Ai £ IP . Let J 73 , C T 3 be the subset of connected horizontal 

cluster forests of scale j containing Aj . Consider the rescaled quantity 

STsS^;^) := nM*^^ J )^( p ))4 ussian (P;G«) 
i=l 

(6.16) 

and: 

-sum over aK F' (£) C {1, . . . , d} \ {*<}, (ij) c{l,...,d}\ 

(JA,r)c{l,...,d}; 

- maximize for (xc), (x' e ), (xa,t)j eac ^ one ranging over its associated 
interval in O 3 . 

CaU ^auTeiian^; A i) the result. Then 

&;a{)<(^t. (6.i7) 

Proof. 

1 . The integral -^G aussian (F- 7 ; G 3 ' 3 ) may be evaluated by using Wick's lemma. 
Each choice of contractions leads, using the numerator in the right- 
hand side of eq. (|6. 12|) . to some term with the correct homogeneity 
factor, X-ov^-^Yl^M-j^NliG^+NliW))^ Consider the resca i e d 

fields d := M 3 ^\. For reasons to be discussed presently, we 
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shall apply Corollary 16.21 to the rescaled fields n(A x ) 7 ^' rescaled (x), 
for some power 7 > 1. 

The possibility to introduce this supplementary scaling factor n(A)~ 7 
comes from the following argument. Split ^ l+M j^ x _ x ^y into 

— , with r = r' + r" , r' , r" > 0. The product 



(l+Mi\x-x'\) r {l+M3\x-x'\Y 

of propagators rifeWFJ) x £i 4> x 'i) contributes, see denominator 

in the right-hand side of eq. (16.12p . a convergence factor K^, ^ ■ 

ElAeFJ ]Ta'~a {d 3 (A, a 0) r ^ 1 ^ or some constant K T » depending only 
on r" . Since the number of intervals A' G IP such that d 3 (A, A') < 
n(A)/4 is < n(A)/2, this means that at least half of the intervals 
A' ~ A are at a d 3 -distance > n(A)/4 from A, so that [jfl 



J] («P(A > A'))~ r " /2 < (n(A)/4)- r "/ 2 ' n ""' = ]p( A )n(A) 



n(A)/2 



-K'r"n(A) 



(6.18) 

On the other hand, taking into account the n(Ai)^ , resp. n(A^) 7 fac- 
tors separated from the rescaled fields in cluster intervals contributes 
]l AeF , n(A)T Arj ( GJJ - A ) < nAeP n ( A ) 7 ° (n(A)) > a Product of so-called 
/oca/ factorials, which is compensated by the above convergence factor 
as soon as r" is chosen large enough. 

We may now apply Corollarv l6.2l eq. (|6.3p to the rescaled fields, which 
yields, using once again N 3 (G 3,3 ; A) = 0(n(A)), 

rf,rescaled/ F j.^jjN < K T, Ae¥ j ™(A) \ K v (F-? ;G") 
Gaussian \ ' ^ — 

0(n(A)) 

0(n(A')) 



n 

agfj 



l + (n(A))" 7 W A ')) 



l + *'(A,A')) r ' 



(6.19) 



Now the sum Y^A'eW (i+dJ(A A')) r/ conver S es as soon as r' > D. Hence 
each term between square brackets is bounded by a constant. Since 
£ AeFJ n(A) = 2|F'| - 2 = 0(\¥ 3 \), one gets: 

jj.rescaled < ^|F^ | ^K f (A^ ;G J J ) /g 20) 

Gaussian — V ' 



16 In D dimensions, becomes iCr^A) 1 /" and eq. (|6.18p holds, with different con- 

stants. 
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2. Associate to a connected forest F J G J- 3 _• and a monomial G h3 as in 

A{ 

(|6.13p its Wick expansion, represented as a sum over a set of con- 
necting pairings of W as in Corollary 16.31 (1), except that iV(A) < 
d(n(A) + ta) + n(A) - the number of fields and half-propagators in 
the M-adic interval A - depends on ¥ J and G 3 ' 3 , and is unbounded 
since n(A) may be arbitrarily large. Hence (to get a finite bound for 
our sum of Gaussian integrals) we shall use the F- 7 -dependent rescaling 
by n(A)~ 7 = C^iV^A) -7 ) of the fields defined in 1., at the price of 
the extension of the exploration procedure described in the proof of 
Corollary 16.31 (2). Note however that the mapping (¥ 3 ,G 3,3 ) i— > con- 
necting pairing is not one-to-one, since a link of the resulting pairing 
¥ may come either from the links of ¥ 3 or from the pairings of G 3 ' 3 ; 
this contributes at most a factor 2 per pairing, hence at most 2 dv / 2 . 

Now, the factor £ A , eD , £?i'=i K((A,i), (A',*')) of Corollary El (2) 
associated to the rescaled fields defined in 1. is bounded up to a 
constant by X^A'eBJ (i+<#(a A')) r ^ °°> nence the result. 

□ 

The above arguments extend easily to single scale Mayer trees of poly- 
mers of scale j. The new rules are: 

(i) there may be some undetermined number of copies of each interval 
A 3 , each with a different color; 

(ii) fields in intervals with different colors are uncorrelated; 

(iii) each cluster forest of a given color is connected; one of them (the red 
one, say) contains a fixed interval, A^; 

(iv) the different cluster forests are connected by Mayer links. These define 
a tree structure on the set of colors, and imply for each link between 2 
colors, say, red and blue, an overlap between one red interval and one 
blue interval (chosen at random if they have several overlaps). 

The proof of Lemma 16.41 (2) is the same as before, except that the explo- 
ration procedure must now take into account Mayer links. Let UMayerO^') 
be the coordination number of a (red, say) polymer P' in the Mayer tree. 
The overlap constraint between P' and its neighbours Pi,... ,P nMa er (P')-i 
in the tree splits into multiple overlaps of order rtj = ni, . . . ,n c > 1 between 
an interval Aj in P' and rij intervals in Ti{ neighbouring trees, IP^ij • • • , IP^^, 
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with m + . . . + n c = UMayeri^') — 1- The exploration procedure at the red 
interval Aj adds Hi to K^, see eq. (|6.7p . corresponding to the number of 
possible choices of neighbouring trees, but Cayley's theorem, see proof of 
Proposition 13.121 yields a factor j- Li ^, < - , 1 - , . Summing over all 

\ n Mayer\^ ) *-)■ n l n c- 

possible values of rtj leads to replacing by Yl ni >o ~1TT~^ = 

The whole procedure must be slightly amended to take into account 
rooted Mayer trees (see §3.4) connecting possibly an interval A G A eX f(P), 
where P G ■p.?-*' i s a polymer with > N ex t, max external legs, to intervals 
without external legs of polymers of type 1. Then one should associate some 
small power of A, say A K with k <C kq (at the price of reducing slightly kq 
in eq. (|6.17p ) to each interval with a field lying in it, while the intervals A 
of the above type and containing moreover no field define intervals of a new 
type (of type 2, say), with no small factor attached to them. The whole 
discussion is very similar to the Remark after Proposition 13.121 For such 
intervals, Ka = 1 must be replaced with 1 + Y^ n >l ^TT! = 1 + 0(\ K ). As 
explained at the end of this paragraph, such a factor is not a problem. 

We may now give a more general, multiscale bound which takes into 
account secondary fields and high-momentum fields. We rescale the low- 
momentum fields by reference to their dropping scale j, and not to their 
production scale k, see §3.3, which leaves outside a supplementary spring 
factor that will be used to fix the horizontal motion of the polymers as in 
§4.1. As mentioned in the Remarks following Corollary 12.71 and Definition 
13.91 we shall not split low-momentum fields 3 into a sum (field aver- 
age)-!- (secondary field) if Pi < —D/2. The following definition is valid in all 
cases: 

Definition 6.5 (spring factors) Let Wipf := d^ipf if Pi > —D/2, with <P 

D 

2 - 



defined by means of a wavelet admitting [Pi + vanishing moments, and 



Pift ■= if ^ < -D/2, so that, by Corollary RT 



^.rescaled^ a )gJ ,rescaled^ ^ } < R (6 . 21) 

(l + d h (A J x ,A J xl )) r 
where ^i> rescaled ^ := M^Ptp^ is the rescaled field, and 

fa = Pi (Pi < -D/2), fa = ft-Lft+^+lJ G [-1-D/2, -D/2) (ft > -D/2). 

(6.22) 
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Example. The cr-field in the (eft, d(j>, cr)-model has ft- = —2a > —1/2. 
Thus low-momentum fields are severed from their averages, and ftj = /3 a — 
1 = -1 - 2a. 

Hypothesis 6.6 (high- momentum fields) Assume 

(i) either that all ft, i = 1 , . . . ,d are < EE' 

(ii) or, more generally, that there is a scale constraint on Ci n t of the fol- 
lowing form: rewriting Ci n t as 

£int(x) = J2 £ £ K^^{x)...^(x), (6.23) 

q>2 l<h,...,iq<dji<...<j q 

then 

( K t;:£ * °) (a, < o, A, + ft 2 < o, . . . , /3 n + . . . + p lq < o) . 

(6.24) 

This condition on the scales of the low-momentum fields is of course 
equivalent to a condition on the scales of the high-momentum fields due to 
the homogeneity of the vertices. 

Note that Hypothesis H6.24H holds true for our (4>, d(j), cr)-model, since 
splitting a vertex leads to one low-momentum field, either d<p or a, with 
respective scaling exponents a — 1, —2a < 0. In general, it has the following 
obvious consequence. 

Lemma 6.7 Assume Hypothesis \6.2J$ holds, and fix I q = (ii, . . . ,i q ). 
Choose e > such that 

e < min (\/3 h \, \/3 h + ft 2 |, . . . , \(3 h + . . . + ftj) (6.25) 
whenever there exists j\ < . . . < j q such that K^'"''- q ^ 0, and let 

li-=--Pi, i = h,...,i q ; 7/:=y^7i- (6-26) 

iel q 

Then ft + 7j > for all i, and 7$ , + . . . + 7j 9 < D for every q' = 1, . . . , q. 
17 which is the case of </> 4 -theory for D > 2 for instance 
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Proof. Since f3 h + . . . + (3 lq = -D, 

7y + • • • + H q < e + (D + fa + . . . + /3 V _J < D. (6.27) 

□ 

Under the above Hypothesis, one has a multiscale generalization of Lemma 
16.41 by considering the contribution of all fields of some fixed scale j. We 
adopt the following convention. If a vertex tjjj (x) is split into i fields of scale 
j and \I q \ — i fields of scale 7^ j, then it contributes a (fractional number) 
of vertices, tj-t, of scale j. This implies a small factor X K ° V for the Gaussian 
bounds at scale j, where v (the total number of vertices at scale j) is a 
fraction with bounded denominator. 

Lemma 6.8 (multiscale generalization) Assume ipi, ...,^C {(9 n V'i)n>0) 
are derivatives of some independent multiscale Gaussian fields ipi, . . . ,tpj. 
Fix some constant kq E (0, 1) as in the previous lemma, as well as some 
reference scale j min < j. 

1. For each k > j, consider a horizontal cluster forest ¥ k E J- k and asso- 
ciated cluster points X(k,x' lh , and choose subsets (Igk), (IL), (^A fe ,r)' ( x A fe ,r! 
as in the previous lemma. Do the same for each h < j, and choose 
for each point X Xgh j X pfa OT X T Ob restriction interval A 3 = 

A{ h ,Ai,A?. h such that A J C A £ h,A' ih or A h respectively. 

Such as choice defines uniquely a monomial G 3 ' 3 as before, and one 
more monomial per different scale, 



£ k £L(F k ) 



n n^vu^* 

A^gF* t=1 

(6.28) 



for k > j, and, see eq. 



for notations, 



Qj,h ._ ^K v(F3;GJ> h ) 



n e e 

e h £L(F h ) AicA (h A'J'CA' 



M- 7l ^- h) Res h Aj 4jx Eh )M V' h) Res h A/j ^ r Jx' £h ) ■ 

(6.29) 



II II E M-^ U - h) ResUi AhT (^,r) 



A h gF h t=1 A^'cA' 1 
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for h < j, where the A J , A are restriction intervals. 

Let v(¥, G>) := v(P; G^)+Ek>j G i ' k ) + 12h<j ^ G jA ) be the 
total number of vertices. Let finally := G^^ Ylk>j G^ ,k Y\h<j G-'^ 
and 

4aussian(^^):= jd^) \\ C$(t* , X g j ; £ , x' &j ) &. (6.30) 
Then 

d 

(6.31) 

where a stands either for a low-momentum scale h < j or a high- 
momentum scale k > j, and uj max < is an in Definition ^. 4\ 

2. Fix the total number of vertices, v := v(¥, GP). Define F^ min ^ {^ min ^) 
as in Corollarv \6.3\ (3). Consider, similarly to the previous lemma, the 
rescaled quantity 

r j,rescaled/ w . ^js 

1 Gaussian ) u ) ■ — 
d 

TT M jft(^(G IJ )+^(P)) TT M^i^ll (¥■ G j ) 

XI 11 Gaussian V. ' r 

(6.32) 

and: 

- sum over all ¥ G J'.w^^™" - *'), (I £k ) C {1, . . . ,d}\{i £k }, (I' £k ) C 
{1, . . .,d}\{i' £k }, (I A fe )T ) C {1, ... ,d} (k > j), and similarly for h < j; 

- sum over all possibly choices of the restriction intervals Aj; 

- maximize for (x f k), (x' ik ), (x^h T ), each one over its associated inter- 
val in D fc , k > j, and for (xgh), (x' £h ), {x&h T ), h < j, each one ranging 
over its associated restriction interval in W (and not B> !). 

Call L^Zf^iv; A jmin ~*) the result. Then 

j^rescaled^. ^j m i n -+j < M -(D+\w* max \)#^m in ^ ^XX K °) V . (6.33) 

Proof. 
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Consider first the factor M -( D +\ u 'max\)#^ : ' mm ^ _ comes f r0 m the part 
of the rescaling spring factors used for fixing horizontally the polymers (see 
§4.1). Let P~ c - *' be one of the connected components of the multi-scale 
forest at scale j, and Ac C A Jmin ^ n W be its intervals of scale j with 
external legs. Then the rescaling spring factors of the corresponding low- 
momentum fields {Til) inext )-+ti- l \x next ), n ext = 1, 2, . . . , N ext (®%~*), yield 

a factor Y^n^l^i ^ M^ n ^t when going down from scale j to scale j — 1, 



where ^ = (3 Hext or p next - 1, see §4.1, are such that £££=i ; ^ < 

D I * 

I ^max I • 

Let us now prove 2. directly, since 1. is a weaker form of 2. Note that by 
the Cauchy-Schwarz type inequality J\fgh\ < (J |/| 3 J \g\ 3 j j/il 3 ) 1 / 3 , one 
may separate low-momentum fields from high-momentum fields and from 
the fields produced at scale j, to which the previous lemma applies p. 

Let us first consider low-momentum fields. We use the same rescaling as 
in the proof of the previous lemma, namely, we consider the rescaled fields 
n(A^)~ 7 5 fc ' rescaled V'i , w ith ^ > 3- The A-?"™" "^-connecting pairing associated 
to (F, Gi) has links of two types: 

(i) links due to pairings (ipl^,), or, more or less equivalently, cluster links 

C\ {in , x e , ; i' £j , x' tj ) of scale j ; 

(ii) cluster links to the propagators (i^a , X£a • i' £a , x' £a ) of scale a ^ j, 
a = k > j in the specific case of low-momentum fields. 

Cluster links of scale k > j (or h < j) contribute a factor ( 1+d fc(A i fc,(A')* ; )) r ' 
which is required both for the bound on P and for that on I k . Since r is 
arbitrary, one chooses it large enough and splits the above factor among 
the different scales of the vertices. On the other hand, the scaling of the 
propagators or is left for the computation of I k or I . Note the 
possible existence of chains of propagators of scale k connecting two vertices 
with low- momentum fields of scale j; summing over all possible chains yields 
the same factor of order ^ 1+d k^l (A') k )) r as soon as r > D. 

By Definition 16.51 the term between square brackets in eq. (16.6P is 
bounded up to a constant (see proof of Corollary 16.31 (3)) by A c i uster (A k ) + 



18 Of course, fields produced at scale j may also be treated on an equal foot with high- 
momentum fields for instance. 
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A low {A k ), where 

4wa») :=EE £ (rnpw at? 

i/=ife'=j A'eD ft ' v v ' " 

(6.34) 

where A k is the unique interval of scale k' such that A k D A fc , and 

d k k' 

A l ™(A k ):= Y YM^'-ny M~^"-j) y 1 , 

v > Z^ Z^ Z^ Z^ (1 + & (A k , A")Y 

i t i'=lk'=j k"=j A"GB fc " 

(6.35) 

The above sum, £ Afc " eDfe » (1+dJ(A ^ Afc » ))r , is of order M D ( k "~i\ which 
yields 

d p k' 

A low (A k ) <KY A£° s , A l r~ R ■= Y M^ k '~ j) Y M&' +D W- j \ 

v ' ~ ^ Pi,Pii PUPii ^— ' < 

i,i'=l k'=j k"=j 

(6.36) 

while clearly A c i uster (A k ) is finite since ft, ft' < 0. 

Let us finish the proof with the assumption that D = 1. There are 3 
different cases: 

- either ft/ < -1; then ft/ = ft/ and Y%>=j M^'+W'-ti = 0(1), 
Efc/=j Af^(fc'-i) = 0(1) since ft < 0; 

- or —1 < ft/ < — |, resp. < ft/ < A: then ft/ = ft/, resp. 
ft/ - 1, and EfeL i M^'+ 1 )( & "^') = 0((A/ - j)M^' +1 )( fe '-i)), Ek#' - 
j)M(ft+ft' +1 )( fe '-j) = O(l) since ft, ft/ < 

- or — i < ft/ < 0, resp. ^ < ft/ < 1; then ft/ = ft/ — 1, resp. ft/ — 2, 

and J2k"=j M^' +I )( fc "-J) = 0(1), while the sum over k! converges as in the 
first case. 

The simpler case D > 2 is left to the reader (there are only 2 subcases: 
/3 < —D/2, or (3 > —D/2, the latter subcase to be split according to the 
value of L/3+ f J). 

Consider now high-momentum fields, produced at a scale h < j (or 
h < j). By the same method, one ends up with the following quantity 



to bound instead of eq. (16.351) : 

d h ti . 

^ > ■ Z^ Z^ Z^ Z^ (l + (ii(Ai,A'J)) r ' 

i,i'=l h'=j m i n h"=j min A'jgDi 

(6.37) 
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where A J , A range over restriction intervals, plus the finite term A c i uster (A h ) 
due to cluster links of scale h as before. Now ^a'j'gidw 

and ^ + jijfii' + 7i' > by Lemma [6771 hence ^^(A- 7 ) is bounded by a 
constant. 

□ 

Note that a small part of the rescaling spring factor M^ k ~^ (ft < 0) 
for low-momentum fields may be used to obtain a factor M~ £ < 1 for e > 
small enough per interval A belonging to a fixed polymer P - in particular 
in empty intervals where there is no field. This simple remark is essential 
in the sequel since various estimates yield a factor 1 + 0(A K ) per interval, 
which may be compensated by this (not so) small factor M~ £ . On the other 
hand, to each vertex or equivalently to each non-empty interval - or to each 
field - is associated a factor of order A K , which may be arbitrarily small. 
This is a general principle for the bounds to come now. 



6.1.3 Gaussian bounds for polymers 

It remains to be seen how these Gaussian estimates, valid for each scale, com- 
bine to give estimates for the (Mayer-extended) polymer evaluation functions 
defined in subsection 3.4. Note that (rescaled) low-momentum field averages 
have been left out; the domination estimates in subsection 6.2 prove that it 
is possible in the case of the (<p, dcfi, <r)-model to bound them while leaving 
a small factor per field, A KO+K o 5 kq,k' > 0, negligible with respect to A K °. 

We take this into account for our next Theorem by choosing a small 
enough exponent kq. 

Theorem 6.1 Fix some reference scale j m i n > and some exponent k > 0. 

LetV J min ^(A j "" n ) be the set of vacuum (Mayer- extended) polymers down 
to scale jmin containing some fixed interval A Jmin of scale j m in ■ Let Eval Ko (P) , 

sum over all multi-scale splitting of vertices into 
cluster forests F J extending over P and monomials G 3 of the product inte- 

grals\[ p j=jmm j d//(^)n te L(FJ) ^fe'^;^'^)^' all fields in G J being 
accompanied as before with the factor A K ° . Then 



E 



Eval w+K j(P) < (K\"y. (6.38) 



pg pJmin (Airmn); ^{vertices of F=v} 

In particular, for A small enough, 

Eval„„ + „.(P) <oo. (6.39) 



E 
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Proof. 

Let •Pj? m!ri_> '(A- : ' m ™ -> ') be the set of vacuum (Mayer-extended) polymers 
down to scale j m in with fixed set of intervals A- 7 ""™ -5- as in the previous 
Lemma. Note that the horizontal fixing scaling factor M~( D+ \ u '™ ax ^# AJmm ^ 
makes it possible to sum over all inclusion links of the polymers. Namely, 
each inclusion link A J C A J_1 - implying necessarily some i-derivative in A J 
- produces a factor M D due to the choice of A J among the M D intervals A € 
ro j such that A C A-?" 1 , which is compensated by the factor M~^ D ^ u ^ iax " 
attached to A J . 



Hence it suffices to prove eq. (|6.38|) for P ranging in the set V§ 



■Jmin 



(A* 1 



Let M 



(A* 1 



Pe"Po mm (AJm<n^);#{vertices of F}=v 



Eval K0+K / ) ( 



Split the small factor per vertex into A K °A re o, and set apart A K ° to get a global 
homogeneity factor X K ° V . Then 



Pe7>, 



o 



E 

" ( A Jrn.in. 



EvaL/ 



(6.40) 



where the total number of vertices is now unrestricted. 



interval A J of scale i. tbr (x,m J) or nl) 



Let P G "p^" 1 " 1- ^ (A jfmi ™ _> ). Consider some product of fields of type q pro- 



duced in some interval A- 7 of scale j, ^,(x AJ t ) or ip Iq \{ it }(xi), interpreted 
as some pairing of tpj q (xf) with ib^x'p), and: 

- choose some non-empty subset I hl 9 h = [i q , ; . . . ,i q ) c I q ; 

- choose some high-momentum scales (ki) ieI hi g h, with j < fc^, < . . . < 
h as in Hypothesis 16. 6\ and restriction intervals (A kx ) i&I high, A* 1 * G ID) fci , 
A 4 fel c A- 7 '; 

- letting I Zo ™ := 

^ < i- 

Then 



I q \I hl9h , choose some low-momentum scales {hi) ieI i ov 



E EEE 


II Res^ 




n ^ 


/Wsh c / A fe, 


i£jhigh 




iei low 



(6.41) 



is the decomposition of ipi into all possible splittings. Any given splitting 
of tpj q is supported on an M-adic interval n ig /h l9 hA fci of size bounded by 



M~ ki i D 



= M-i D ■ M~ {k w~ j)D . . . Af-^-Vi) 13 

< M" iD M" (7 V + '" +7i ' )( ^ ri) . . . M~ liq<ykiq ~ kx i-^ 

= M~ jD M'^^i hl9hll{kl ~ j ' ) (6.42) 
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by Lemma 167 

Fixing ki, letting j range over all scales < ki and changing notations 

ih,j) ->■ (j, h) yields the spring factors M~""e h ^~ h \M~^ I 't h 3 ~ h ,M~ 7/ aV (j ~' i) 
of eq. (|6.29|) . The remaining factor M^' D = |A-?| may be rewritten as 



M 



■3 Hi 



which is distributed between the different fields, if)', 



-0: 



j'i,rescaled 



M 3 '"^' as in §6.1.2. Recall from the end of the preceding paragraph that 
each interval of P comes with a factor Ad £ <C 1. Hence, letting jmax 

be the 

maximal scale of a given polymer P one has by Lemma [ 



( A 3 min 



EvaL 



< 



p 



Jmax 

^ n 

Jmax — Jmin J — Jmin 



M~ 



j,rescaled/ »j rai „- 
Gaussian V > 



«>1 



Jmax — Jmin 



if A if chosen small enough so that in particular (with the constant K of eq. 

§33}) E t ,>i(^ K °)" < K ' XK '° < 1 - M_e - 



□ 



6.1.4 Combinatorial factors 

The last point - for the Gaussian part of the final bounds - is to control the 
combinatorial factors due to the horizontal and vertical cluster expansions. 
Let us show briefly how to do this. 

As a general rule, differentiating a product of n fields yields n terms 
(by the Leibniz formula). This implies supplementary combinatorial fac- 
tors when estimating the polymer evaluation functions -F(P), compared 
to the estimates of Theorem 16.11 Consider the 0(n(A J )) derivations in 
a given interval A- 7 due to horizontal/ vertical cluster expansion at scale 
j. A field produced by one such derivation may be acted upon by an- 
other one in the same interval, yielding a local factorial (0(n(A jf )))°( n ( A: ')) . 
Otherwise, derivations in A J act on fields produced at an earlier stage 
in some interval A k C A J belonging to the same polymer P as A J . In- 
tegrating over these, and borrowing some small power of A from one of 
the differentiated fields, yields at most an averaged factor in A J of or- 
der K := |^7|A K y^ fc>i y^ AfeePnDfc . Afc(=AJ f Ak dx. Once again, there are 

0(n(A J )) such derivations. One may always gain a local factorial n ^jy. 
to some arbitrary power (see proof of Lemma [63]); using ^,L, < e K , and 
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multiplying over all scales j and all intervals A J £ P n IP , one obtains 

expA K ^ J2 M ~ ik ~ j) ' ( 6 - 44 ) 

k A fc gPnB fe j<k 

hence a factor of order 1 + 0(A K ) per interval A £ P, compensated by some 
factor M~ £ as explained at the end of §6.1.2. 

6.2 Domination bounds 

Unlike Gaussian bounds, which are rather sophisticated, these are essentially 
based on the simple fact that |x|e _j4 l x l = A^^A^e'^ < KA^ 1 if A > 0. 

Lemma 6.9 (domination) Let ip be a multiscale Gaussian field with scal- 
ing dimension (3. Then 

\(Ti/))^ k - 1 \A k )\ n exp-\ K M mf)k ■ j^J (iT4>)^ k -^y {x) dx 

< K n n n/ mX -Kn/m M -npk_ ^ 

Proof. 

Let u := (T^)^( fe - 1 )(A fe ) and v := ^ f Ah ((T%f>)^ k -^) m (s) dx; by 
Holder's inequality, |u| < v 1 /" 1 , so that 

n \ n n J 

< K n n n l m \- Kn l m M- nfik . (6.46) 

□ 

Example ((0, c?<^>, cr)-model). Lemma 16.91 implies in particular the fol- 
lowing four kinds of low-momentum field domination: 

(i) Av£ 4<< 5£ 4 terms 

Consider low-momentum fields a produced at a scale k by letting some 
derivation 4- or dt (due resp. to horizontal and vertical cluster expan- 
sions of scale k) act on £ 4 . When p — k is large enough, they will 
be dominated by the part of the counterterm 8C4 which is coupled to 
bP- 1 . Eq. (RT351) yields (using 1 - t 2 > (1 - t) 2 for t £ [0, 1]) 



(l-t k A r 



ni 



n e -A 2 (l-(t!) 2 )M( 1 - 4a ) fe / A fc \(Ta)^( k - 1 )\ 2 (x)dx 



< (KM 2ka n-t ) . (6.47) 
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Note the factor (1 — t k A ) n coming from the rest term in the Taylor- 
Lagrange expansion as in §3.3, with n = N e xt,max + 0(n(A)). This is 
precisely the reason why we chose to Taylor expand to order N ext: max + 
0(n(A)) and not simply to order N ex tmax- The other terms in the 
Taylor-Lagrange expansion have t\ = 0. 

Replacing in the exponential \ 2 M^- 4a ^ k by the term ss A 2 M( 1-4o! ^ _1 ) , 
one gains a supplementary "petit facteur" (^M~2 (i-^Xp-l-*)^ _ 

Av£ 4<(5£l2 -terms 

When /) — A; is too small, the "petit facteur" is not small any more. 
One dominates by some fraction (say, one tenth) of the boundary term, 
jq5Cu instead. Again, eq. (|Q5j) yields (using 1 - (t k A ) 6 > (1 - ^) 6 ) 



(l-t k A ) n 



\{Ta)^ k - l \A k ) 



H e -^A 3 (l-(ti) 6 )Md- 12 «)«= f Ak \(Ta)^-V\S(x)dx 



< 



[K\^M 2ak Y < (K'X^M 2ak y 



Replacing in the exponential M^ 1 12a ) k by the term M^ 1 l2a )f> present 
in 5C\2 (note that 1 — 12a < !), one loses this time a supplementary 

(12a — 1 / 7 \ \ ^ 
M 6 yp~ k ) J . However, it is accompanied by a small 

factor Aa per field. 

Assume p — k = 0, 1, . . . , q. We fix q such that X^M ql2 °e = A 1 / 4 , i.e. 
1 = 2(i2a-i) lT \nli^ • Thus the "petit facteur" in the preceding Av£ 4< $c 4 
case - where p — k > q - is at most A K , where k = ^i 2 ~a-i) • 

Av,5 £l2<( 5 £l2 -terms 

Consider low-momentum fields a produced from some vertex A p by 
letting some derivation act on 8C\2 this time. It produces i < 5 low- 
momentum ex-fields, accompanied by A 3 . Again, eq. (|6.45|) yields, by 
dominating by one tenth of the boundary term, j^dCu, and leaving 
out once and for all the i-coefficients, 



(6.48) 
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1 

n! 



X 3 {(Ta)^ {p '- 1) (AP')) 



< (K'X 3,i - 1/2 M 2a P' ' M( 1 *Tr 1 H p -' / ) 



(6.49) 



The corresponding vertex produces after integration a factor \A P ' \ .M ( 12ct = 

M- l2a P.MP-p\ multiplied by M 2a f>' = M 2a PM~^ p - p ^ M' 122 ^^-^ 

per high-momentum field, and by M 2ap> M^ir 1 ^-^ X 3 ^ 1 / 2 = M 2ap M~^- pn > X 3 ^' 1 / 2 

per low- momentum field. Thus all together one has obtained a small 

factor < A 1 / 2 per vertex, to be shared between the six fields, and 

M e ( p ~P ) < 1 per high-momentum field. 

(iv) Av 8 c 4 <8c 4 

Consider low-momentum fields a produced in an interval A k of scale 
k by letting some derivation act on dC^ or on the counterterm of scale 
p, — ^(tx) 2 {(Ta)^ p {x)) 2 . It produces at most one low-momentum 
ex-field, accompanied by b p ' = (If /X) ■ X instead of A. 

Again, eq. f|6.45f) yields, by dominating as in (i) by the part of 
the counterterm SC^ which is coupled to b p ~ 1 , a factor 0(M 2ka ■ 
M -I(i-4aO( P -fc)) per field, alas multiplied by ¥>'/X « AAf'* 1 - 4 "'. The 
rest of the argument goes as in subsection 6.1.4 - a general argument 
called " aplatissement du fortement connexe" in (colloquial) French. 
The factor M~^ 1_4a ^ p_fc ^ may be simply bounded by 1. 

Such terms may be produced only in intervals A p ' C A k such that 
A p £ P. Integrating over all such intervals - and taking into account 
the M _2fca -scaling of the high-momentum field a k left behind - yields 
at most XK := X Y^ P >>k M~ 4a ( p '~ fc )#{P D W'} per low-momentum 

field produced, all together (A 1 / 2 )™ • (X^K) n . One may always gain a 
local factorial using ^j- < e K and multiplying over all scales k and 
all intervals A fc £ P n B fc , one gets 

X v ' 2 • exp A^ ]T ^ ^ M-^p'-^ (6.50) 

where v is the number of vertices where such low-momentum fields 
have been produced, hence a factor Aa per vertex and a factor of 
order 1 + O(Aa) per interval A £ P. 
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6.3 Final bounds 

The main Theorem is the following. 

Theorem 6.2 1. There exists a constant, positive two-by-two matrix K 
such that 

V = K\ 2 M^- Aa \l + O(A)). (6.51) 

2. There exists a constant K' such that the Mayer bound of Proposition 
\3.12\ for the scale j free energy, 

\f^ p (X)\<K'X(l + 0(X)) (6.52) 

holds uniformly in j . 
Proof. 

Let us first prove eq. (|6.5ip for b 3 . 

Consider a product (G-monomial) x (product of propagators) as in sub- 
section 3.3, written generically as GC. Multiply it by a product of aver- 
aged low-momentum fields of one of the four types Av£ 4<< 5£ 4 , Av£ 4<< 5£ 12 , 
Av5£ 12< 5£ 12 or Av,5£ 4<( 5£ 4 , see subsection 6.2, generically written as Avi ow . 

We make the following induction hypothesis: 

Induction hypothesis, b = K(l + 0(A)) for all k > j for some scale- 
independent constant K, where b k := A~ 2 M -{ l - ia ) k k j s fo e rescaled mass 
counterterm of scale k. 

We shall soon see how to compute the constant K. For the time being, 
we must bound 

j dn a {<t>)dn a {a) E Mo W GCe- S \n^ +5c ^ sc ^ < ^ dx • 

FeV^(Af) (G,C) 

(6.53) 

where |P| is the support of the polymer P with two external <r-legs, in the 
notation of Proposition 13.111 

First step (domination of the mass counterterm). 

Note first that the term between square brackets in eq. f|6.53j) is negative 
when X := \M~ 2pa \\{Ta)^ p \\ is small. Up to unessential coefficients, it is 
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equal to M P (A~ 3 A" 6 - {t p x ) 2 X 2 ), which is minimal, of order M p \ 9 / 2 - for 
tx ~ 1, which is the worst case - when X is of order A 3//4 . The factor t x 
in front of (Tc)^^ -1 ) selects the intervals in D p belonging to the polymer. 
Hence the exponential in eq. (|633]) is bounded by e K MPxV2dx ^ a \\ 
together a factor of order 1 + 0(A 9 / 2 ) per interval A 6 PnB p of the polymer 
with t p A / 0. 

Second step (domination of low-momentum fields). Split Av[ ow e I\p\( Ci+SCi+ 2 SCl2 ^ • 
from the expression in eq. (|6.53p . As shown in subsection 6.2, these pro- 
duce a small factor of order A K , k > per field. More precisely, any 
k < inf (j, ^2a-i) ) * s suitable. There remains (by using the Cauchy- 



Schwarz inequality) to bound XV g -p^(aj) S(g,C) (/ dfj, s (<f))dfx a (a)\GC\ 



2\h 



Third step (computation of V s ). 

Let us now estimate V by means of our induction hypothesis and of the 
Gaussian bounds of subsection 6.1. Consider for instance the diagonal term 
lr, , . The terms with the fewest number of vertices are: 

- the term with vertex obtained by applying twice to the coun- 
terterm 5 £4, namely, lP + , (by definition); 



Figure 5: Main part of the mass counterterm of scale j. 
the polymers with two vertices, see Fig. \5\ which sum up to 

A 'E / C^(x,y)(-d 2 )Cl i (x,y)dy =: X 2 M^~^K V , 



k>j 
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where 



-+, VKoo k ■.= m-^ 1 - 4 -) | ier 4 v (opcode = / i£r 4 V(0x^ 2 (£K, 

(6.54) 

a scale-independent quantity. As in section 1, the non-diagonal counterterm 
iP + _ or b J _ , may be computed in the same way, yielding in the end a scale- 
independent positive two-by-two matrix. 

More complicated polymers are of order APt 1 '^ \ 2 g(\ 2 V , X 2 (b k ) k>j , A), 
where V := A -2 M - -^ 1-40 ^ is the rescaled mass counterterm as in the 
induction hypothesis. The function g is C°° in a neighbourhood of and 
vanishes at 0. Hence, by the implicit function theorem, V = KM^ 1 ~ 4:a \l + 
0(A)) as in Proposition 13.111 

The bound for the scale j free energy f 3 ^ p is now straightforward. □ 

Bounds for n-point functions are easy generalizations of the preceding 
Theorem. Consider for instance the 2-point function (| J r 4>i(£,)\ 2 }\, with 
Mi < |£| < M J+1 . By momentum conservation, and by definition of the 
Fourier partition of unity, see subsection 2.1, this is equal to the sum over 
ji, j2 = j,j i 1 of (F(j} 3 i {^T^i {— £))a- The term of order in A is given 
by the Gaussian evaluation E J-^(^)J-<^(— £) . Further terms involve 
at least one cr-propagator with a small factor (see Lemma 12. 13j) of order 
inf(l, A/J( fe 1 ; 4 " ) ) < K'm£(l, A~ 2 M~( p_J ^ 1 ~ 4a )) which goes to when p -»• oo. 
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